۳۱ ۲ ۲ د‎ AUR AD EM I-A NAUK 
INESE Y ۰۳۴ 1 12 ۲ KT 


E ACTA PHYSICA 
POLONICA 


۳۵۰۰ Pn DU 


DWUMIESIECZNIK 


Orders and inquires concerning 
Acta Physica Polonica 


— complete sets, volumes and single fascicules — 
as well as other 
Polish scientific periodicals 
l published ۱ E 
before and after the war, d 
regularly and irregularly, ۱ 
are to be sent to: - 


Export and Import Enterprise ,RUCH" 


Warszawa 1, P.O. Box 154, Poland | X 
Ask for cataloques, folders and sample copies. cow 


Vol. XV (1956) ACTA PHYSICA POLONICA Fasc. 4 


EXPANSION OF SINGULAR FUNCTIONS ASSOCIATED WITH 
THE KLEIN-GORDON EQUATION 


By L. Inretp and J. PLEBANSKI 


Institute of Physics of the Polish Academy of Sciences and Institute of Physics 
of Warsaw University 


(Received March 1, 1956) 


Two alternative approximation methods involving expansions in powers of c^! 
or c are given for the Klein-Gordon equation. By using these methods, expansions of the 


more important singular functions such as 4, A, and A( are obtained. The connection 
is shown between this procedure and two alternative methods for formally integrating 
the Klein-Gordon equation on the basis of analogies with the Yukawa:and the harmonic 
oscillator equations. The main purpose of this paper is to prepare the mathematical basis 
for the subsequent application of the Einstein-Infeld-Hoffman method to the مد موہ‎ 
of motion in meson ۰ 


§ 1. Introduction 


Functions of the type A, AO, A, A p» etc., (e. g., see Schwinger 1949) associated 
with the Klein-Gordon equation occupy an important place in the mathematical 
formalism of modern theory. Their particular importance is connected with the fact 
that in both classical and quantum theory, the proper choice of Green functions and 
propagators is determined by the physical character of the interaction. 

On the other hand, the new approximation method (Einstein, Infeld, Hoffman 
1938) is well known from general relativity theory. Roughly speaking, this method is 
based on the expansion of all quantities in the theory into a series of expressions of 
different orders in c^! (c is the speed of light). Simultaneously, the orders of the initial 
expressions in the expansions are postulated for each of the particular physical quanti- 
ties. This turned out to be the only effective method for deriving the equations of 

‘motion for particles interacting gravitationally from the field equations (e. g., see Ein- 
stein and Infeld 1949, Infeld 1954). For a long time, however, it has been known that 
this method permits the successful construction of Green functions in problems in- 
volving the d'Alembert equation in classical electrodynamics. This made possible the 
simple derivation of the equations of motion in electrodynamics (Infeld and Wallace 
1940). 

(207) 
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Recently, Werle (Werle 1952, 1953,) applied this method to the Klein-Gordon 
equation. He obtained interesting results, particularly on the level of quantum theory 
(Werle 1953,, 1953). 

The purpose of this paper is to throw some new light on the applicability of the 
new approximation method to the Klein-Gordon equation, and as far as possible, to 
explain the mathematical basis for its application. In order to do this, we will develop 
two approximation methods, one connected with the expansion in powers of l/c 


mc : é : : 
where u = a: is regarded as being of the first order in l/c; the other, connected with 


a series expansion in powers of c. This allows us to obtain the expansions of all the 
more important singular functions and to discuss their significance. In § 5, we com- 
pare these expansions with the results derived by certain formal methods of integrat- 
ing the Klein-Gordon equation. Expressions similar to our symbolic formulae in § 5 
were used, for example, by Nambu (Nambu 1950). 

The expansions obtained in this paper are not merely another heuristical path 
connected with an approximation method leading to the singular functions of the 
Klein-Gordon equation. We intend to use these expansions in the equations of motion 
of classical meson theory (Infeld and Plebański). We also believe, however, that these 
expansions may be useful wherever linear electrodynamics and mesondynamics 
touche on non-linear general relativity theory. So far, the only effective method 
permitting analysis within the framework of non-linear general relativity theory is 
precisely the new approximation method. 


§ 2. Definition of Singular Functions Associated with the Klein-Gordon Equation, 
and Their More Important Properties 


In this section, we will discuss the definition and properties of the singular func- 
tions which will be needed for the treatment in this paper. 
We are interested in the relativistic, invariant, real solutions of the equation: 


(E - 2?) G (x?) = — by (X) (2.1) 

The solution inside the future and past cones for the point «*=0 (regions T+ and 

T^) can depend only on the invariant quantity c = (x, د(٤8‎ = قم)‎ t? — r2)*. Out- 
d 


side the cones (Region S) the solution can depend only on the invariant parameter 
o —(—x, (1۶ = (r? — گا (ئوقم‎ 


df 
Inside, T*, T- and S, function G must satisfy the equation!: 
Evite 1 § 1 
for T+ and, T-; E پت ال‎ + (12 - i) 1G (2.2a) 
di^ 9 0 0 1 | 
ES (+ (| «e -o (2.2b) 


1 We express لا‎ in terms of T or c. 
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Because of (2.2), it is possible to represent G inside the individual regions as: 


inside T+: G = at 1 4 P (س)‎ (2.33) 
pt 11 
^ inside T-: 6 — a-2! ae TU. Dum (2.3b) 
inside S: ES (2.39 
m po 


where a*t, ره‎ d’, bt, b-, b are arbitrary constants. Indeed, these equations-do not 
tell us anything about the behaviour of G on the cone. We have written them merely 
to see clearly where and how Bessel functions appear in the singular functions of the 
Klein-Gordon equation. 

The basic singular functions are defined as follows: 


f ی‎ a 2,315050 
نت سے‎ na 
J; x*>0 
Z ole) ۀ‎ (9) + e E " art sz cis S A 
۱ IPSA 2,232 2» 0 
40 — EL. TT (2.4c) 


= ب‎ e) pa x0 


Table Y contains the more important integral representations of tee functions 
` together with their definitions. 
| These functions satisfy the equations: » 
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We use the 4 to define the function 40۰ 


Ao = Ag (& — Zt, t; xis. cL 252/128) 
0 


It can easily be shown with the help of (2.5), (2.6), (2.7) that: 


a موی‎ 5 
رو سےا‎ Am = - کا یذ‎ A)  —0—4 40| _ (2.9) 
t=to 
Equations (2.5), (2.7), and (2.9) permit us to write explicitly the diete of the 
Cauchy problem for the equation (O — p?) u (x, t) = — 4ro (x, t) in the form: 


+~ 
u (x, t) = 47+ | d x fa Ao (x — x’, t, tst) o (x, t) + (2.10) 
1 , 9 = =, , AX =y Qu وم‎ 
- 2 | ae EE (t — to, x — x) u (x, to) + At -2) 9 gag |. 


This equation very adeptly expresses the meaning of A and A, respectively, 
as a Green function and propagator of the Klein-Gordon! equation 
Solutions of the Klein-Gordon equation representable as an advanced or 


retarded potential have particular physical لا دو‎ The fundamental ralations 


are: SEC 


A => (Art 4 Ando, Lt gad — Are (2.11a) 
so that: 

1 + سم و ord‏ = 

Aret = A — 7 A, A“ = A+ pi 4 (2.11b) 


We will call all solutions of the equation ((J—y?) u = - وہ‎ which have the 
form: 


u (x, رز(‎ = dove 1 dyx' E dt’ At, (t—t, x—x) E) (2.12) 
retarded (advanced) potentials. This terminology is suggested by the fact that: 


T کڈ - )4 وي‎ F') سے‎ .)2.13( 
Ee Att Ir ۳ سو و(‎ 


7 ۷ 0: t<t’ 


Are (£ — t',7 e — i 2 races 


! The covariant generalization of (2. 10) for dealing with the initial و‎ on any T o does 


not present any difficulty. 
` ` **) We will consistently use a dash over all solutions of the inhomogeneous equation: 


(D — y’) G = ô, (2°). 
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Aa ES سے‎ 4i 2 راچ سیر‎ 
ado => , حتف‎ LAN. 9 .ےی گس‎ - 
Aado (t — t'r —7') = ie 28 AEK t',r—r') (2.13b) 
_a 0; t 7 ti 1 
لا‎ 4 (t — و غو( کہ ہہ‎ 
which means that only those values of the sources o (r’, t’) for times t’ earlier (later) 
than ¢ enter into (2.12). 
If we now let u go to zero in all the equations written so far, we obtain the relations 


for the basic singular functions of the d'Alembert equation. We will tabulate the more 
important ones. We define 


D= — Lie = م د‎ mgl. | د )+ +( و‎ (2.14a) 
af 4n 877 cr c 
ارک‎ ART 0 enu je a 
و = )6.27 ۵ (د) مرچ و‎ 8 (e r?) nes 0 darte lt 2 
(2.14b) 
pius ہے ہت‎ ee eee 
q xat 2 ct کم‎ Ancr PUT ARA [2.14c) 
c c 
The representations of these functions are given in Table II- 
These functions satisfy the equations: 
OD = -èq (x), OD=0, OD = 0. (2.15) 
D and D are related by: 
= 1 olet : 
For function D the following holds: 
92 1 ۵ 
D = U = — MES =s : 
1-0 ) 31270| تی سب 5 ای‎ 


The function: 


n n - l[t,—t' T 
Do) = D Ex tus bo) تہ ہے 085 نن‎ x’ : 
(0) ) (x — x 0 : و"‎ ae i=") 96 - t'x—x) (2.18) 


has the properties 


- 


| 8 = 
Ux Do) = — ô, (x* — 9, ظ‎ | = 0 = = Doj _ (2.19) 
=f 


ot 


Using Dass we can write àn equation analogous to (2. 10) as a solution of the 
` Cauchy problem for the d'Alembert equation. 
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Functions D“ and D are related to D and D by equations analogous to (2.11): 
D= - عق‎ LDN رر‎ D 7 (2.20a) 


A = = = 1 
D“ =D — 2D D” — D+ x (2.20b) 
D“ and D^? are the kernels of the retarded and advanced potentials of the 
d'Alembert equation. 


83. Application of Approximation Methods to the d'Alembert Equation 


Our actual problem is to examine the singular functions associated with the 
Klein-Gordon equation from the point of view of two alternative approximation methods 
(involving expansions in powers of l/c or c). For greater clarity, we will, however, 
first consider an analogous problem in the much simpler case of the d'Alembert 
equation: j 
1 ۶ 
E وع‎ 

The considerations in this section will be elementary; they will not give any new 
results. They will merely serve to illustrate the structure of the two approximation 
methods. 


Ou (x°) = — dao (x°) ; O = v? — (3.1) 


We will look for the spherically symmetric solutions of the equation: 


1.3? la 
[v - i2 G = — کو‎ ۵ (0. (3 2a) 
This equation can obviously be written: 
o? > 
c? y? — E G = — c ô (x) ô (t). (3.2b) 


G can be expanded into a power series in negative or positive powers of c: 


1 CES V 1 ۱ 
G-A-— وه‎ +2 a ES (3.3a) 


6-8-۰ dete) b am, عو‎ 
2 1 


1=0 1-0 


Inserting (3.3a) into (3.2a) and requiring that this equation be satisfied for every 
order in 1/c, we obtain: 


Via = پی۷۷٦‎ va = 0 
ae (3,4). 
a; Lol. 


^ (3.4a") 
رو‎ Ot? و-رو‎ s a 208 47 
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Similarly, inserting (3.3b) into (3.2b) and requiring that this equation be satisfied 
for every order in c, we obtain: 


92 9 

یہ ۱ وۀ  <‏ ہے 

A (3.4b’) 22 (3.4b") 
کبس و‎ b ll . سے‎ = DLE 

Ot? و رو رو‎ 9t? on 44 21-1 


We will call the approximation methods connected with equations (3.4a) and 
(3.4b) 4 method and B method, respectively. 

Equations (3.4a) and (3.4b) were deduced from (3.2a, b) by expanding in powers 
of 1/c and c. These equations however, could have been obtained without taking c 
explicitly but by using x? = ct and assuming that for A method, the derivatives 
with respect to x? are of smaller order than the derivatives with respect to x*; for B 
method, the opposite is the case. 

We will first consider A method, and, in particular, equations (3.4a’). The equation 


1 
for the zero order is satisfied by a = — ô (t)!. 
o 41 


(QI 2 
We puta = 220 f(r); for f we have: f = ٠ ; g ات‎ f This equation 
2 47 2 2 0 T r 0r ولهو‎ 


= 
is satisfied by se . Thus, the set of functions 
Jer ۱ 
راو سل‎ po) 


will be a solution of (3.4a^). In particular, (3.4a") will obviously be satisfied if we take 


a = 0. This with the help of (3.3a) and (3.5), we have obtained a solution of (3.2a) 
21+1 
exclusively in terms of the even a, namely: 
1 


ES OUT 


Assuming, however, that we can use a کے‎ s expansion for the ۵ function, we have: 


3s Ln xo 2040. Le که‎ eee) 8 
شر وس‎ i! (e Y 1 "A E 


1-0 


= aos [oe 2) 6-2) 


1 It serves no purpose to add toa a solution of the homogeneous equation v? ^ — 0 restricted by 


0 
the assumption of spherical symmetry to any function of time. We will start with ath a solution in equa- 


tions (3.4a”). Treating the function of time as part of a would disturb the breaking up of the solution 


into powers of ۰ » 
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Comparing this with (2.14a) we therefore have: 
A” = D. (3.7) 


The above can be regarded as a simple heuristic deirvation of D by the obtaining 


of its expansion. What is the meaning of this expansion? In practice, D, as a Green 
function, always appears in expressions of the type: 


Rp — 
U (3, t) = 4nc | dt 1 dy DOLE ۵ x): (3.82) 


Introducing into the above equation the expansion (3.6), changing the order of the 
summations and integration, and making use of the properties of the ۵ function, we 
see that 


eroe Jere- 0000۷۷" 


i-0 


Certainly, it makes sense to speak of the convergence of this series only when the 
integrals in the individual expressions exist. This, indeed, occurs when the sources o 
are either singularities of the type à, (x— £ (t)), or vanish rapidly enough, e. g., expo- 
nentially, at infinity. 

The case of singular sources is particularly important in problems concerning 
the derivation of the equations of motion. When the integrals in the series (3.8b) 
exist, the question of the convergence of the series can be raised. The formal formulae 
(3.7) (i. e., formal from the viewpoint of ordinary analysis, which does not recognize 
the ۵ function), show that if this series converges, it will converge to the value of the 
integral in (3.83). 

We may add that (3.8b) could be obtained directly as a solution of (3.1) without 
using D, by applying A method to (3.1), i. e., by taking its expansion as a series in 
even powers of l/c. 

We now pass to (3.4a"). The spherically symmetric solution of Y? a — 0 is an 


۰ c 2 > 1 
arbitrary function of time a = YW (t). Beginning from this function, it can easily be 
1 ۳ 


seen that solutions of equation (3.4a") will be 


pl 9? ۱ : 


The series 41444 constructed from these odd a 
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1 € 1 i [2l 2۱ 1 
بر — په‎ B ا سے‎ | 
8 : 2 ریزو‎ 09773 ٥6 2 11-17) ij v (0° ہج‎ )8.10( 


is a solution of the homogeneous equation [] 4444 = Ol, 

If in the case of 14, the initial term was almost uniquely defined by the condition 
of spherical symmetry, then the arbitrary function y (t) would appear in 4٦44. The 
question arises: With which function y (t) are we to begin the approximation procedure 
for A°“4? The answer is simple: we begin with the function which we would like to 


obtain for r = 0:lim godd = . Y (t). 


f0 


Thus if we want to obtain 4944 = D, then we must choose 


Aer "۷ئ۷ ےا‎ 8 ۷ 
p= lim د ۴ رر‎ aC ۳ On 


If we want to obtain 49% = D™, we must choose 


y = hm c? — neo" = 7ن0‎ 
ae peo 22 (coL ete) 2772 t2 


It is easy to see that this is so. When y = x 6M t), then by (3.10): 
T, 


1 v put 1 E 
27 cr La (QE + 1( 77 Pe 


1 1 
BET 010) |) = E = gez 1۵ 6 rà — ۵ (t — rjc)) 


D = A, (3.11) 


Aodd — 


and 2 (2.14b) :? 


1 [f in integrating equation (3.4a”) step by step, we were to add, after each step, a solution of the 


homogeneous equation, e. g.. 9 (t), then in -quation (3.10), instead of y (t), there would appear 9 (t) = 
214-1 


9 (t). By not including these Q (t), we not do suffer from 3 loss of generality. 
THT oy 21+ 1 


1 
ip^ 
9 
x 


2 Having the expansions of D and D, we immediately obtain the expansions of Dret and Dade; 
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Starting with y = — ET “ we obtain from (3.10): 
I ee AONE lode tr ان‎ 
ہے - - مھ‎ Y سے‎ Fa تر‎ 1c 0721 23 EI 
2z?c £e (20 + 1)! ۱۵۸ (9E. co 2720 fe t c 
1 1 


22 cet? — r?’ 


: 
where the last sign of equality is valid only when ^, و‎ < l, i. e., inside the cone. 


2 
۲ ۰ 
Outside the cone, i. e., when—> Ap > 1 the series 4°% has no meaning. Thus, method 4 


gives us: 
2 
DW = Aod for 72 1ج‎ (3,12) 
We will now take up method B, and, in ELT equations (3.4.b") Since 
P 5 | ۶ | = ۵ (t), the solution for : will be ; = — ۳ ô; (x). This suggests that we 


look for 6 in the form de = g (0) Es ô; (x) where g= * ||. Substituting this into 
21 


(3.4b^, we find that © me 
2) 21-2 


These equations are satisfied by g — 54 BITS : 
Thus, b will be given by:! 
2 
کی سے‎ E و‎ 


Taking b = 0 as a solution of 7 4b”), and using (3.3b), we obtain as the solution of 
21+1 


(3.2b) the series composed exclusively of even b, namely: 
1 


1=0 


mw چو‎ t GG (ati 
C= Be تار‎ OIM dep" &, G) (3.14) 


Representing ô, (x) in the form 258 ili dskei* * and, after changing the order of summa- 
oo 


2 2 2۲ 
1 By و8 لاپ‎ (x), we understand simply the result obtained by operating with ( Z s + ج‎ T z) 
on the function ۵ (x) ô (y) ô (z). 
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tion and integration, and reducing the series to an integral function under the integral 
sign, we obtain: 


1 1 1 >> sin kct 
Be هه‎ aa و‎ ikx 
2117 ix e k 


~ 


Comparing this with Table II, we see that 


D = B”. (3.15) 


Returning to (3.4b"), we see that the equation ET b — 0 is satisfied by e = 
= t f (r) + g (r) where f (r) and g (r) are arbitrary. مات‎ from this : we can SEIS 
see that the remaining equations of (3.4b’’) are satisfied by the fauc ae 


2+1 9ع‎ 
FEDORAS (3.16) 


2041 
It is obvious that the series B°% composed entirely of odd b in the form 
1 


~ 


21 سے پچ 
f(r) + c? (ct) v?! g (r) (3.17)‏ لپ Eu. BaF E Bi‏ 
£i un ; QUE I y: De: |‏ 


is a solution of the sical equation c? ح‎ 4 — 0. 


If the initial term of the series B” was given almost uniquely, B°” should then 


contain the two arbitrary functions f (r) and g (r). The question arises of now to choose 
odd 


: oar 0 
these functions. The answer is simple: since 4ھ‎ = ¢ g (r), and 3; = 
۶-0 t=0 


c? f(r), then f and g are determined by the initial conditions for the function we wish 
to obtain. For example, if we wish to obtain 944م‎ = D, then me to (2.17), we 


should take g = 0, f= E (x). Then, 2 TS TT A? à, cE 


If we take 6, (x) in a Fourier representation and change the order of integration 
and summation, this ۶9۶ then reduces to: 


ME TE sin ket- 
B و‎ ey] bet 7 


Comparing this E Table II. 2, we have: | 
D وت‎ (3.18) 
į 


In order to obtain B°% = D, we should choose, according to (2.14c) 
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We then have: 


TOES ہے‎ 


r2 +2 ر2‎ 2 c?;2 — r? 


22 1 
The last sign of equality is valid only when e > 1, i. ہم‎ outside the cone. 
T 


Method B therefore gives us: 


2۶2 
DO — Be for 2 <1. (3.19) 
T 


Inside the cone, the series is divergent. 

We will now sum up the above discussion. Method A, which is connected with 
the expansion in powers 1/c, and method B, which is connected with the expansion in 
powers of c, both work in conjunction with the singular functions of the d'Alembert 
equation in the following manner: Methods A and B can both be used to construct the 
series expansions of D and D in the entire space-time continuum. In the first method 
however, it is necessary to assume that the Taylor formula can be applied to the 6 
function, while in.the second, it is necessary to use the Fourier representation to show 
that the expansions are convergent. Thus, in both methods, the individual terms of the 
series for D and D are constructed from the 6 function, so that the convergence criterion 

of ordinary analysis cannot be applied. As to D®, method A permits its expansion 


Table III 


Method A Method B 


1 
t b 
D) زود‎ pers 8 1 1 (ct)?! ۱ 1 | =0 8 
1 c= و اج‎ amit ۳ 
و نار ۶ > : سب‎ 
7+1 : re کہا‎ 
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inside the cone only, while method B, outside the cone only. These expansions are made 
with ordinary functions and are convergent in the usual sense. 
The expansions obtained are tabulated in table III. 


84. Approximation Method, for the Klein-Gordon Equation 


In this section, we will generalize the procedure of 8 3 for the case of the Klein- 
-Gordon equation: 


(O - p?) U (x°) = - Ame جصو)‎ = 7 . (4.1) 
We must, therefore, develop the approximation procedure in a manner somewhat 
parallel to method A und B for the d'Alembert equation. In this connection, we come 
to the following question: What order of l/c should be taken for u in method A? 
It would seem the simplest to assume u to be of the order zero!. 
This assumption, however, leads to mathematical difficulties connected both with 
the convergence of the respective expansions, and the possibility of constructing the 
functions A, A, and A“ from these expansions. 


In this paper, we will use a variant of method A in which u is regarded as a magni- 


a 


tude of the first order in l/c, i. e., u = q/c, where q = is regarded as quantity 


independent of c. Why do we make this assumption? First of all, this method, as we 
will soon show, allows the procedure of § 3 to be generalized in a simple and natural 


way, and permits the construction of the expansions of A, A, and A“. This genera- 


0 
lization will be based on replacement of the operator 5 of §3 by the operator 


23 
E 
view in 85. Another argument can be taken from the general theory of relativity. 
In using the Einstein-Infeld-Hoffman method to derive the equations of motion, we 
must assume that the mass is of the second order in 3 
Since we intend to use in the general theory of relativity the expansions obtained 


+ q?. Some new light will be thrown on this matter from another point of 


cm 
here, it is therefore logical to assume the second order for the m of ہر‎ =— as well 


h 
as for other masses appearing in the theory. 
In using method B for the Klein-Gordon equation, it turns out that m should 
be assumed to be of the zero order in c. The discussion in 8 5 will cast some new light 


1 This assumption was made by J. Werle in his papers [loc. cit.] 
? Newton's equations have the form of mass x second derivative of position with respect to time = _ 
mass x mass x a function of position. In the E. I. H. method, the time derivative is of the second order 
in l/c; therefore, in order to use Newton's equations in the fourth order, the mass must be taken to be 
of the second order in 1/c. 
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on this matter. Thus, in the end, we will seek spherically symmetric solutions of the 


equation. 
(DI — 45) G (x^) = — 9 (9) (4.2) 
which can be written in the equivalent forms 
۱ EE > ۵ (t) 
(for method A) |» 2i E 4- 2] G = — و۵‎ (x) یسوی‎ (4.3a) 
92 
(for method B) E (v? — u?) — A G = — cô, (x) ô (t), (4.3b) 
as a power series in l/c or c: 
ا‎ EAS 1 
ود 277 یيو‎ (4.4a) 
G= B =e NA bel te b. cH, (4.45) 
Tn اس‎ 17 


Inserting (4.4a) into (4.3a) and requiring that (4.3a) be satisfied independently for 


each power of l/c, we are led to the equations: 


via = — à, © ô )( vio = 
0 1 
32 , 9? 22 
Dem RE 2 4.5a یں‎ sd Pact 21 a 4.5a 
x 7 E Tid ۳ ) ya (3 9 ۳ ( 
1>1 LO 


Similarly, inserting (4.4b) into (4.3b) and requiring that (4.3b) be satisfied independen- 


tly for each power of c, we arrive at the equations: 


7 (4.5b’) 8 (4.5b”) 
و سب‎ =(y?2 — y2) b مخ ات سام‎ ٢ 9 7 
23 2l (v ۸ و‎ Qt? ie (v ٨ ۳9 

SI‏ لتد 


We will first discuss method A, and, in particular, equations (4.5a’). It is obvious 
1 
that 2 rpm ۵ (t) satisfies the equation of the zero order. Employing the same line 


of argument which we used in 8 3, we can easily show that (4.5a’) is satisfied by the set 
of functions determined by the zero step approximation: 


i 1 m [1‏ 92[ ابا 
om (te) 0= > (ijesme ao‏ و 
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The following question now arises: What is represented by the series 4”, which 


contains only these even ے‎ and which is a solution of the inhomogeneous equation, 


namely, the series: 1 


Sie Ly, 1 1 
و ور‎ [DNA 271 D HE ie ari 9) 


In Appendix I it is shown that 


0 ۶۶۶ء0 
At? = 0 cou 4.8‏ 
rart > 0 rD‏ ی جع + )0( 
This function, which satisfies the equation ( O — u?) 4% = — ô p(x") does not resemble‏ 


any of the normally considered singular functions of the Klein-Gordon equation. It- 


was obtained by a direct generalization of the procedure used for D in the case of the 
d'Alembert equation. We introduce the following notation for this function: 


۸۳ = 4 (4.9) 
df 


Function 4۰ is equal to zero inside the cone for the point x^ = 0 and differs from zero 
outside this cone. As shown in Appendix I, this function i is connected with 4 by the 


relation: 
APA ٢ 5 ۹ )4.10( 
where 
سک‎ 3 (PD AEE Se, (4.11) 
— Ar \ L (uo) | uo; xax" <0 


` Asa result of the appearance of A?” (which, as can be easily seen, is a solution of the 
homogeneous equation ([]— u’) 4:4- 0) in the expression for 1, a reciprocity 
. arises between properties of 4” and A. Indeed, when u — 0, then 4144 — so that 


both ۸۳ and A correspond to D. From (4.10) we conclude that the function 4 can 
be expanded in powers of 1/6 only when we know the expansion for ۸۳, This expan- 
sion will be obtained by using the approximation method in relation to the solutions 


of the homogeneous equation, i. e., (4.5a"). 


. From (4.5a' ^), we see that as a result of the condition of spherical symmetry, we can 
put a = vy (t) where y (t) is an arbitrary function. As in $3, it is easily decuded that 
1 


the solution for the subsequent steps will be the functions: 


ral 92 سو‎ 
uc (21 4-1)! [5 9 PO od )4.12( 
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The solution of the homogeneous equation obtained by the approximation method will 


thus have the form: 


io 1 SAE. ا ا-‎ 2 Ls (4.13) 
"9۹۹ ۳ ei U ۱۷٣٣ ; 


Here, as in § 3, the answer to the question of how to choose y (t) is as follows: Bearing 


1 ۱ : 
in mind that lim 4747 = — y (t), the function y (t) is determined, for 0, by the 
c 
r=0 
value if the function which we want to obtain for 4°“. In particular, we are interested 


in the three functions; و4۳49‎ A and A“. The corresponding v (t) are easily calculated: 


lq 

i MA‏ 7۸0 ۵ تہ 

y () = lim ctdett — پر کے کے‎ (go (4.14) 

)4.14( ال د RC ETE‏ که کا 

1 ae |t] 7 و‎ 27 

: ۱ G 

= 2 (1) —— Hn 1 

y )۵( = lim è AM = و‎ FN, (qo) (4.14) 


We will now show that if we substitute ۷ (¢) from. (4.14) into (4.13), we obtain 
for 1944 the functions ۸۶4۵, A, and A“. In proving this, we will pay particular attention 
to the manner in which the respective series converge. 


We will first take y (t) from. (4.14’). On the basis of (A.8e) of Appendix II, we 


have: 


)( = کے سر پګ‎ + | 00 4.13 
Therefore, using (A.8d), we have, by (4.13): 


r 21 1 22 l 
PEE (iat) vo- 


1=0 
1 نو‎ ۶ 1 Pp 141 
= c ایل ما تیل نتم‎ c 2 hr 
472 3 1 )27 +- 1(1 "e ed Jo (qt) — (4.16) 
aI 
~ 2aci 7 2 (3) (2) Ji+1 (qt). 
1-0 
Comparing the last series with the Lommel expansion (Watson 1944, p. 142), 
رد‎ co (—)™ (4? LU 1)” 1 m 
RU Ed pees 35] r+ (2) (4.17) 
m=0 


1 The respective limits can be calculated either b 
and A) (i. e., (4.11), (2.4b), and (2.4c)), 
and Table 1.4). 


y using the explicit representations of 1044, A, 
or on the basis of their Fourier time-representation (i. e., A: 3. 
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which is valid for any arbitrary complex 4 and z, we find that 4ء‎ converges, for every 
t and r, to the quantity 


1 RE رام‎ 
{oda سا لجا‎ 2,2 - à 
Arc yer- rà Ura J (alc Vc t T ). 
Returning to the relation u = g/e and distinguishing between the cases x,x*2 0, 
we thus obtain: 


| 


۰ 1 بت‎ 1 PE 141 
ota = 27, y 7 3 (2) Ji+1 (qt) = 
7-0 
BNF (adler همه‎ > 0 
47 


۱ 


|= Aoda 
۱ lh(uc)uo xax“ > 0 j 
which is as should be expected. 


(4.18) 
If we now make use of the expansion of 1 (4.7) and 41۰44 (4.18), we can write 
the expansion of 1 in powers of l/c with the help of (4.10): 


= 1 
و )0 اد‎ oda ٣ 


9? d 

sat *( ô (1) + (4.19) 
1 سے‎ 1 1 q 1+1 

ge نو‎ z Jr (g). 


This expansion is valid, without restriction, for every r and t. For q > 0, it obviously 


corresponds to the expansion for D. (For q > 0, 1۶44 vanishes). If we put this expansion 
into 


aos سے‎ > — > 
UG dore 1 dt’ f شوه‎ A (t—t', x—x) o (t', x’) 


(4.20a) 
and employ the properties of the 6 function , we obtain, after changing the order of 
summation and integration, 


€ ' 


1 
1 n |x — x71 o (t, x) + 


q 11 : : - 
9 (t = 5 | 1+09 [t "De 0 , x). 

i (4.20b) 
If the integrals in the individual terms of both series exist, and the series converge, 
then U of (4.20b) will equal U of (4.20a). Indeed, the first part of U of (4.20b) can 
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easily be obtained by applying the approximation method directly to the equation 
(O-p?) U = —4no. Asa result of the previous discussion, it is necessary to subtract 
the second series which contains the typical meson ‘tail-in-time’. This must be done 
since we wish to obtain the solution associated with A and not with ۰ 

We will now insert y (t) from (4.14”) into (4.13). We note, first of all, that in accor- 


dance with the relation ic Sod m 6 (t) and (A.8e) we can write: 
0۸ 2 ا‎ 
2-0 0+ نوا‎ IE ده‎ 2 
ie E RLS n u^ Mc TEILT کو‎ t), 4.21 
so that: 


ri 3 ۹ 1 r21 32 à D, " ME 
Be OLSSDI تو‎ 9 y) a 1 og + 9 [el 0 (qt) (4.22) 


For what follows , we will have to evaluate the a given above. This can be done in 
21+1 
two ways: 1) by differentiating directly, and using the formulae of Appendix II and 


the properties of the derivative of the sign function; 2) by using a method of induction 
with respect to 7. This result was obtained by direct calculation, but it is incomparably 
simpler to check it by induction. The result is: 


t 1 F2 ۹ 1+1 
He TT J 41 (gt) + 
1 
1 E I+% 
+ بو‎ E E (4.23) 
k=0 


Thus the corresponding series for 4°% has the explicit form: 
MUSS وم کر‎ aloe p 
Aodd — A ES بخ‎ ۷ 
یی‎ R 8 1 3 Jia (ge) + 


1 w/(r\” 1 ٤ ۱+۷۵ . 
+ One 7 2 ۰ GI M q» ( k ) ۵)2+1-20( (2). (4.24) 
k=0 


We must answer the question about the convergence of this series and its sum. 


The sum of this series is most easily obtained by using the Fourier representation for a; 
namely, starting with the integral 12+1 


و0 
Jo (qt) = 3 sin (g | £ | chy) dg (4.25)‏ 
0 
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by a change of variables, we arrive at the relation: 


t = i - 2 -iot شه‎ OD E 
ا‎ 0 = [+ [am Sem چي‎ E 


Using (4.26) in (4.22), we easily find that! 


(p — qy 9 eie (4.27‏ گے PE (—y | + | do‏ ا و و 
a 70 [oo q E 27)‏ 
D EE‏ 


Constructing A“ from these a, we have: 
27+1 


oda — roy. 
Cc d +1 سے[‎ 


+۵2 کٹ‎ es پا‎ q?) +% e -iot بد‎ 


7 — 2 l 
cM i 
"e. 


7 X ٢ 21+1 
1 NL af jes verom 7 7: Vo? — q rie وو‎ 


Assuming the series converges, we MX the order of the summation and integration. 
Comparing (4.28) with the representation of the function A in (L4), we conclude 


1 By a change of variables, the integral of (4.27) can easily be put in the form: 
co 


Gated t= dg (q sh y)2!+2 sin (qt ch (و‎ 
= 5 sin Ci ۸ 
وا‎ 01 peto 


It might have appeared that by using the well-known integral: 
oo 


2i z/ 2 
0 


1 1 1 
which is valid when 0 > arg z < x and Rep > us qr when arg z — 0 و وم دي‎ aD de ir 


(4.23) could easily have beeri obtained. However, this is not the case. The representation for a is actually 
21+1 


the limiting „forbidden“ case for this integral. According to rigorous treatment of analysis, this inte- 


gral is not defined. This is clear from the presence of the derivative of the 6 function in a. 
21+1 
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that if the series for 4774 is convergent, then: 


T ود‎ Lay شو‎ A eae (4.29) 
47 La تب‎ Quer ae TY] ۳ i : 
This last proof, although simple, has many shortcomings. It does not give any 
information about the convergence of the series for 4°44. In Appendix III, however, 
we present more a complicated proof which contains this information. 


Thus, we have already made use of the expansions for functions ا‎ and A in powers 
of l/c. These expansions were obtained with the help of the approximation method. 
It is obvious from (2.11b) that we can now write the expansions for A’ and ۰ 

There still remains the question of function 4P. In order to obtain this function 

و 1 


with our method, we set, according to )4.14'''(, y (t) = 27 9, à (q t). We note, 
however, that by (A.8e), we can write: 
1 0 1 /| 
y (t) = 25,5; ^ (qt) — 2 E 3 2 No (qt). (4.30) 


Using (A.8d), we then have: 


ye‏ يهو [S AM OE.‏ یں اص 
کے yy rre sed AOE FT ae)‏ رت 


] (aM 
22 Il (2) Nisa (qt. (4.31) 
By (4.13), we form A? from these a. 
20141 
وا‎ 2 1 ۸ 141 
odd. تب‎ m b SEE q 
4 3 2 7 3 2 2 c! a 9t Ni (qt). — (432). 


But for any arbitrary Bessel function, the restricted Lommel expansion is valid (Watson 
1944, p. 143), namely: 


~ 


Z, (Az) = x m CP gn E = ٦ Za a 000122970 a 


m=0 


Comparing (4.32) with (4.33), we see that if r2 < ¢2¢2, then ۹۹١ from (4.32) converges 
to the value: 


È Mla ۷۶ - اھ‎ 
int ya e 


Returning to the invariant notation, and noting that u = q/c, we compare the 
above result with (2.4c). Thus, 
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E ۶۲ٔ۰ء"۶‎ 2 N, (ut) 
Aodd — X سم هی ]ہے‎ PATE nr 
ہت‎ L! (5) B x C Ar ut d SES 
=0 
7200۰۰۳ 

If x, x > 0, i. e., outside the cone, series (4.32) has no meaning. 

We will now discuss method B, i. e., (4.5b),. There arises the question of whether 
this method is capable of giving us the expansions of the functions A, A, and ۰ 

We first consider equations (4.5b’). By generalizing the procedure of § 3, it is 
easy to see that solutions of these equations are: 


ELI 0 4.45 
21 AE - 1)! às (2). (4.45) 
Solutions of the equation (c?[] —c?u?) G = — ۰0 (x) ۵ (t) formed entirely from 
even b will therefore be: 
i 
1 بپ‎ (ee : 
eU سه د‎ bai T 2 

Ms e * Tin بش‎ prse Di K?) ۵ (x) . (4.46) 


We will use the representation Ee E, dak eikX for ۵ (2 After changing the order 


of summation and integration, and reducing the series to an integral function under the 
integral sign, we obtain: 


>> 2 
Be — ; Mp sin VE? + کي‎ Tas 
Vk? + p? 

Comparing this with Table I, we see that: 
Be = ۰ (4,47) 
eem to (4.5b") and generalizing the procedure of 8 3, we see that these equations 
are satisfied by the functions: 

t2 1 


2 2 د‎ 
re DISCI وو‎ Qr ay (^ - wa Q) (4.48) 


where f (r) and g (r) are arbitrary. It is obvious that the series poda composed entirely 
of such odd b, and having the form 
1 


Bodd —c b للم‎ +1 - 


mper E 


1=0 


is a solution " 7 کہ جو‎ equation (c? — c?u?) B^? = 0. Similarly as in § 3, 
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odd 
we note that 64ھ‎ — 2 oS = c? f(r); thus f and g are determined 
t=0 t=0 
by the initial conditions of the function which we wish to obtain for B^? For 
example, if we beg that B” = A, we should assume, according to (2.7), that 


ge f s وت‎ Then: 


[21 


Vt POE 
د‎ (2L + 1)! (v? u’) ۵ (x) c 


Bodd — 


which reduces, after using the Fourier representation for Ó( (x) and changing the 
order of summation and aes to the relation 


Bodd — با‎ faker hk نی‎ zen ct. (4.50) 
ET Ba), Vk + ڈیر‎ 
This is the representation of function A as given in Table II, so that: 
ped — A (4.51) 
If we desire to obtain B^? — A, we must, according to (2.4c), take 
1 Ki) . 


RO E DEE 
From Appendix IV, (A. 26), we have, however: 


0 
Oe r ۔‎ 9g?c? v? ٢ u’) Ko (ur). (4.52) 
Thus, on the basis of (A.26), we obtain for the n of (4.49): 
1 ارت‎ 1 1 
روان لا ہے‎ e y OP >0 
= 1=0 

z (ct)! u 1+1 

=a (E) 77 (4.53) 
=0 


This series can be written, however, with the help of Hanckel functions: 


/ 8۱ےھ 2 = 0 ۱ 1 

edd E E و ا‎ ee: ( Jace iur 

Bas. t = Ly ) 0 2۳ (ur) . (4.54) 
Comparing (4.54) with (4.43), we see that B°% is a restricted Lommel series, which con- 


cr 
verges, if cum e l, to the quantity: 


2,2 
: رن ا‎ e | 
Bodd — LA stg gU eini2 H PD V T? ae 


Expansion of Singular Functions 93] 


Returning to the function K and comparing the above results with (2.4c) we see that 
when x,,x^ > 0, we have: 


1 x^ (ct)! P ۶ کا‎ (uo 
p = 1 Y €) (5) س جح ریت‎ = AD. (455) 


If x°x, < O i. r., inside the cone for the point x" = 0, the series (4.53) has no meaning. 

We now collect our results together. Both method A and method B are applied 
to the singular functions of the Klein-Gordon equation in a parallel manner. Expan- 
sions of ا‎ and 4 in the entire space-time continuum, can be obtained by either method. 
However, the usual criteria of convergency cannot be applied to these expansions, since 
they contain the derivatives of the ۵ function. As regards the function AU, method A 
forms the expansion entirely within the cone, while method B, entirely outside the cone. 
These expansions are composed of ordinary functions and converge (or diverge in the 
respective regions) in the usual sense. 

If we wish to obtain A using method A, then, as has already been stated, we would 
be forced to introduce new relativistic singular functions ۸۳ and 1۹۹4 which would be 
different from zero outside the cone. Whether these functions can find use as kernels 
of the interaction integrals remains an open question. 


Table IV 


Method: A 


1 ۶-1 1 / ۶ e) 

9 سا ا‎ a. ô (t) + 

Anc X 9211 c?! 5 v 0 
PU 


ir 2 1+1 
1 1 r 9 
^ aad NEA a Jisa (9; 
1-0 


Method B 


1 zt = (ct)2/+1 
17 
1-0 


al 


(y? —u 3)1 ا‎ (4) 


1+1 21 ہم 
Am r 1 |‏ 


ud sla‏ 2 اوس 
ô; (x)‏ — — 
ub à‏ سس فو 2 4 
س« شون xd sy e‏ 
k-0‏ 
eet‏ | 20 
Ni t‏ — د 
(qt)‏ 1+1 8 )+ 


#4 < 0 
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We may add that in classical theory, only 4 and A appear, so that the fact that ۳ 
can be expanded into a series in l/c only inside the cone, does not have any signifi- 
cance here, as we are using the approximation method connected with the expansion 
in powers of 1/c. In quantum theory, AO, AF, etg., play an essential role, and, as seen 
from our discussion, it is not so simple to carry out the expansion in powers of l/c, at 
least not in that variant of the approximation method where w is of the first order. 
In other expansion schemes, this case is the most convenient one from the mathematical 
point of view. 

The expansions obtained in this section are collected in the table IV. 


§5. Symbolic Formulae 


The relations arrived at in the previous sections have a characteristic similarity 
in structure. In this section we will show that this similarity results from the fact that 
it is possible to develop two symbolic methods for the singular functions associated 
with the Klein-Gordon equation. 

Before proceeding to more general considerations, we will set down a few characte- 
ristic facts. 


1) According to (A.1) and (A.2), the function ۸ can be written: 


+~ 


1 سوب یرہش و توس[‎ 
ےا میں یج یں‎ —io کے رم‎ 2 ics NL IER —iot — 
A که‎ | dee ۶ ch — - Vg? o EIE ch 7 y^ + 28 2 due =i“ == 


—^—- 


)5.1( تچ تد دو کو 
Be ZEE‏ 00ج "goo oy t+‏ 


This symbolic formula has a definite meaning. The operator ch Vy ۳ a 
x 
can be interpreted in two ways: 
+o, 


: ر ہے‎ Q2 M 
d) as MEE مر 2 اہر‎ em , b) as an operator xb makes f (xo) = [dve-*"* g(r) 


يم -- 


a function with the transform ۷ u?— v? g (v). Since ch z is a series in even 
powers of z, the sign of the radical presents no problem. 


2) According to (4.13), the general solution of the homogeneous equation can 
be written: 


2 


: cx 1 21 1 
odd — ہے‎ 7 . 
i: 2 یف‎ 7+2 21 + 1 (= ad : y (t) — 


1-0 Dg 
Te 


= f dt' y (t) R(t ره‎ «( (53) 


-m~n 
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where 


E 7 1 [982 í 
60-2 QI DI 7 (re 1 903 
120 


2 
1 shr pw? + ES 
a 0 
idem B ۵ (x). (5.3) 
2 دس‎ 
La 22 


This last symbolic formula can either be regarded as a compact notation for the series 


formed by differential operators, or interpreted from the point of view of Fourier 
analysis. In the latter case: 


EL Vez. doc nc 
ELSE TEE xa 
T 


sh- ۷ تو - قو‎ 
T 1 pa d em 
TCT 


w?‏ — رہ 


r +~ 
1 1 ch uw y q? — @2 
zc rane E 
927631 r 74 
0 axe 


Comparison with (5.1) shows that: 


R (t, ( سے‎ 


r 


R (t, r) = — | r’ EL (t, r’) (5.5) 
T$ 
According to (A.21) and (A.19), 
0, DESC | 7 | 
1 t 2 
Ae (t, r) = | — (96) 
8x 7 8r |, 0۱۷۵۵-۲۸۱ 
Therefore $ 
0, ] > | t | 
وا‎ = a 


5. 
h (q Vr/c? — t?) r>c|t| - 
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کے 


Equation (5.2) allows us to write the solution of the homogeneous equation obtained 
by method A in the form: 


ib 
a a +0] تہ‎ AR (5,8) 


The operator pet ee appears in the symbolic formulae (5.1) and (5.3). In real 
Xo 


cases, however, only the square of this operator occurs. 
3) It follows from (4.49) that ۶ can be written symbolically in the form 


8پ e Y C‏ هم هې 


1-0 


Exe n تسا‎ f (1) + cà cos (ct Vu? — v?) g (r). (5.9) 


The operators appearing here can be interpreted in two ways, either as a series 


c )) 4 
Bodd — مه‎ See 
5 بے‎ QI + (1 


of differential operations, or by their action on Fourier transforms. For example, if 


Try fd ket hr g (k) then 
cos (ct V u? — ۷2( fC) = f ds k cos ct V u+ k g (à). 


The operator which is characteristic for method B, namely Vu3—v? appears in these 
formulae. In real cases, only the square of this operator occurs. 

All these operational formulae find their roots in two alternative approaches for 
carrying out the formal integration of the Klein-Gordon equation: 


)5.10( .95( رہ 2 - =€ کم - (O‏ 


These methods are organically connected with our two approximation (expansion) 
methods, methods A nad B. Equation (5.10) can be written formally in two ways: 


2 
E ۸ et | G = — ôs (x) ô (xo) (ه5.11)‎ 


ps qs | C= ۵۰۸ ba. 6.11) 


The form of (5.11a) suggests that this equation could be integrated by the method used 
رن‎ 
2 Ox > 
In doing so, it would be necessary to regard the functional dependence of x° in the 
inhomogeneity and in the arbitrary constants as parameteric. The functions of the 


for solving the inhomogeneous Yukawa static equation with a “mass” w= 
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2 
77 r , 3 
mass u = 2 + 
7 9x? 
functions of x?. 


The form of (5.11b) suggests that it could be integrated by using the method for 
solving the inhomogeneous harmonic oscillator equations with a frequency of œw = 


obtained in this manner would be operators acting on these 


= Vu?—v?. Proceeding in this manner, the functional dependence of x? in the inhomo- 
geneity and in the arbitrary constants, should be regarded as parameteric. The function 


of the frequency’ w = ۷-٥ obtained as the solution would then be operators 
acting on the functions of x. 
It will be more convenient to consider first the method suggested by (5.11b). 


First of all, we should define what we mean by the operator Yu—v?. We will define 
it by its action on the Fourier transforms. Thus: 


F(Vg?—v?) | dye g (k) = f dk eî F V + 12) g(t). 612 


The basic property (Vu?— v2)? = u2—v? will obviously remain in force. 
We will look for solutions of (5.11b) in the form: 


G = cos xy Vu?— V? A (xo, x) + sin xg Yu— v? B (xs, x). ۱ (5.13) 


; d d 
The variation of constants method leads to equations for —— A and ——B: 
dxo dxo 


dB _‏ سم کم کد 
VP — y? TA‏ وه 2 سا پا Vu?‏ مد cos‏ 
Vp? — y? sin x, Vu! — v? = 27 — y? cos x, Vu? — y? — 694 (x)‏ — 
dA JB‏ 
which when solved for ST and —— me , leads to:‏ 
Xo‏ 
اس مہ a‏ 
ax i u? = y?‏ 
dB _ COS Xo Vu* — V^ yu? — ۶‏ 
س (x) = ۵ (xo)‏ و۵ (x9)‏ ۵ 
dx ۷ vi Vu? — ۶‏ 


Integrating these equations, we obtain: 


ihe sg 
A= A GB - S IS 


EO) اند‎ | Vu? 
Thus for any arbitrary 4, (x) ada By (x), we e the solution: 
G= cos zy Vu? — A? A, (x) 4 sin مد‎ Vu? — y? By (x) + 
1 ; visage تی‎ 
a, sin xy Vu? — v 4, (3). 


i el - Vots vt 


à, (۰ 


it rei (x). (5.14) 


(5.15) 
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We now compare this with (5.13) and with the representations in Table I,2 and con- 


sider special cases for وغھ‎ and ۰ 


a) Putting Ay = By = 0 we have 


پور و _2 s‏ 
(5.16a)‏ ترخا وو سے ته اد د توب ماو په بل pe‏ 
|xol Vu? — y?‏ 2 


b) Putting A, = 0, B, = — " 63 (x) in the solution of the homogeneous 


u*—v? 
equation, we have 
در‎ ss 2 
ا هت وس‎ es Sa ار ا‎ ME (5.16b) 
Vu? — v? 


1 ۱ 
c) Putting Ay = Van 84 (x), By = 0 in the solution of the homogeneous 
equation, we have: 
Ve 

E وه انر‎ EAE (5.16c) 

Cy 

۷۶ —v 
Formula (5.16c) requires some comment. If the action of V p?—V? is regarded 


as action on a Fourier transform, then this symbolic identity is valid without restriction. 
However, if according to the discussion at the end of § 4, we regard (5.16c) to be: 


2 > 
4 = M E ny og 0 = 
1-0 
= x xo?! 2 2۱7 4-1 1 
z I! (uris eae — ہے‎ Ko (ur) 


then this relation is valid only outside the world cone of the origin of the coordinate 
system. 
The formal method described above for solving the Klein-Gordon equation is 


particularly convenient in problems involving initial conditions. We note that the equa- 
tion 


di ۱ 
ET 2 U (x) = f (xo) راه‎ (5.17a) 


with the initial conditions مد‎ = Xo —-U-—U (xg). dU. du (x) has the solution 


ےتا 


1 t 
1 Tt can easily be shown that سح‎ 6, (2) = E Ko (ur) 
22$ 


(Vu? — y?) 
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Xo 


v= Xo. مت‎ 00 + cos v (xo — xo) U (xo) ^R 


Xo 


sin v (xy — x) dU — 


4 : يد‎ (5.17b) 


Generalizing the above, we see that the equation 


02 ہے سے تفص‎ 7 ds ۱ 
E a + (p? - zd U (xy, x) = oto (xg, x) (5.18a) 
0 
with the initial conditions 
— 
U| =U a, = TU s, s) 
AD Xo | وت رد‎ ^ dx 


has ihe solution 


Xo 


> i 2 2 يي‎ m d > 
U (o, 2) = 4a f AE 
` و‎ LV 


وھ 


p لع وت۔۳ هه سس در‎ Gd). (5.18b) 


Using (5.16), we easily verify that (5.18b) is equivalent to (2.10). 
From the above derivations, it is clear that the same expansions obtained by 
method B in the previous section result from the above symbolic method when the 


integral functions of the operator Vu3— v? are expanded into a power series. 
We will now consider the symbolic method for solving the KJein-Gordon equa- 
tion suggested by (5.11a) and the analogy with the Yukawa ES equation. First of 


all, we should define what we mean by the operator u? ugs — If we require 


E 
merely that y” سا‎ mr 
Xo 


defined by the relation: 


: 9? hd -inx 
Vet ite Ve gal m ° g (v) = 
8 (J +f ) dye E (v) Vy? — p? g (») + (5.19) 


٢ 


+ | dy e» € (v) Vu? — v? g (v) 


77 


, the action of this operator can be 
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In (5.19), ۷۷ رم - 2ہ‎ 2 and Vu— v? are understood to be the positive roots and é(») is 
an arbitrary function which can take on only the values —1 and +1. 


2 
If we require that lim] / 42 + m l 2 we obtain the hint that it will be 
اق‎ 020 


می 
the interval (—oo, — u) and (u, oo). As regards‏ وک = (ر) ۶ convenient to choose‏ 
y‏ 
the interval (— u, + u), we do not have any hints as to the choice of £ (v). It turns‏ 
2 
out that it will be convenient to use for pet T both operators. In this case,‏ 
20 


é (») will have a constant value of —1 or +1 in the interval (—y, + y). 
Therefore, we define two operators: 


(+) AS. سے‎ A 
y^ + جح‎ | aem g (v) 7 —i f+ Í dy ae e-9* g (y) + 
۲ ۸ 
=u 
(=> +~ tu ~ 
y^ 4 ۳ سعما‎ g»)-2-i i 1 dy i V»? — u? eir (v) + 
9x? df 2 y: 1 
ات ^ بت‎ te 
— 1 dy Vu? — ڈو‎ e-irs g .)م(‎ ue 
- 


Both operators have the required properties: 


(-) ۹ )+( 
۶ 92 ?9 ا 02 
سا = ان سے 2 
بت )7 ھن 0 — مدق ین y- 1 zm‏ 


: 9? 9 8? 
li 2 - ed: 2 
1 u? + SES Tusc lim] / u? + 273 (5.21b) 


Using these operators, equation (5.11 a) can be taken in two alternative ways: 


(*) 
د‎ ۰ 23 59 Vs 
V م‎ ۳ aua] | & = — 56 d(x). (3.22) 
But we know that the equation 
(v5 ۸۵ *( 6 == ad, (x) (5.23 a) 


has the general spherically symmetric solution: 


l ck u'r hw 
ILS SM sh ur 
وسو دد اد ووس‎ (5.23b) 


1 
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where ۵ is an arbitrary constant. Comparison with (5.23a) shows that equation (5.22) 
has a general spherically symmetric solution in the form: 


(+) (+) 
1 2 1 / 8? 
سد ہے‎ 20 '- 09 2 EE. کت‎ 
6 chr u? + TE ۵ (xy) + : shr u? + !يچ‎ 0 
z} (5.24) 
l 32 . 
a رت‎ 2 E. 
F-shr]/ we axe, 8 Mo) 


where f(xy) and g (xo) are arbitrary. Since ch z is even, it does not make a difference 
which sign is taken for the radical in the first term. Comparing (5.24) with (5.1) we` 


identify the first term with ۰ 


(ae) 


- 1 1 
Ae = 7 ch TER ٢ تو‎ 0 (oon = T= ch y» 2 "UMEN و‎ Ô (xo) = (5.25) 


(>) 


Wine SE CET 
Haar د‎ ALLAY 32 Oa 


In order to interpret the remaining terms, we keep in mind that f and g can be 
chosen so that that part of G connected with the solution of the homogeneous equa- 
tion can be interpreted as other singular functions. 


Thus, from (4.16), we have: 


1 = وسوا‎ 1 3 1+1 
odd — یا یت یی‎ 
á ET Ed Jo (qt) = 


=% mas me + و‎ 3h + (2.26) 
1 UM NU en 
r Vn ات‎ Vato EC OEE 


Similarly, because of (4.22): 


1 v 9 pu 1 د‎ = 


r E ید‎ ^A E " 
Aster yip, T+ $a ۷ P+ aq) 9 c20 
l ۷ 24 9* 7 24 9* X, 
موي‎ LH u 3 2 u 57 2 | | Jo (H o). 
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In these formulae, it makes no difference which sign is taken for the radicals since 
they appear only in even powers. We note, however, that: 


Ye. So. ۵ (xo) = 
32 
2 EX 
Hs 9x? 
ioe in 
روم‎ e- P» 
= سا‎ -— dv dol dv = 
3 if یی رغ سه‎ fT 2 ese 
|| 
-~ a -u 
~ 1 u 
5 1 dy مد‎ Rl 1 10ت رز‎ )5.28( 


y? — gy? Vu? — »? 


i 


1 
پچ‎ fa sin (ux, chy) + a IE cos (ux, sin g) = 
0 


0 
: 
From (5.26), (5.27), and ae we have: 


5 (4 + مسر‎ ir t Ea P TA vd و‎ T Jo (4 xo) = 
2 امد۱‎ 
1 9 1 
E 2 : 
gr t a "ct u ye ی‎ E 
23 )5.29( 
wer s 
20 


1 M 9? 
rag ase H^ sca 09) 


Adding and am cs these two equations, we have: 


cs ۵ (5.30a)‏ - هم un "yn‏ + 860ج بر + هم ۱٢‏ و 


اوت مان سا Gh)‏ 


1 
odd — 0 


Comparing this result with (5.24), we see that the functions A and A? are obtained 


from (5.24) as solutions of the homogeneous equation by putting f — L ۵ (xy) 


Expansion of Singular Functions 24] 
= ARSE = مد‎ Th i 
& = و ہے ٥ا جر‎ (<o), g اا‎ xz (xo). The other symbolic formulae are 


easily btained by making use of (5.25) and (5.30). For example, since E = 


1 
رو‎ 21044, we have: 


— 1 1 
وب‎ (cr ID my En ۵ (xo) + 
1 1 تر‎ 
ET E : E hr u? + L س‎ eT. mu ô ) (xg) = (5.31) 


wl 2 is 2 
ہے او رو‎ L? +H تش‎ Ta? + exp u? ٢ 27 ۵ (xy). 
Further, since aks Sie z A = Av — b 7 + -) we will have: 
: )+( 
— 1 
بس ور‎ VE 2 due 
Arer Tema La mpg تت‎ s shr 7 / u + £s ۵ (xy) 


5 + و‎ 7 ME = laa) 


(+) 
El exp( = 6 zn) 3) (5.322) 
Ager Ox? : ۱ 
اوځ‎ 
1040 — به‎ ep + u? + pu ۵ (xo) (5.32b) 
Azer 9x? : | 


Because of (5.21b), equations (5.32) obviously correspond to the expressions 
for Daw when مر‎ > 0. 

It is easy to obtain the operational formulae for the function A“. To do this, 
we use the representation of this function from Table I.4. We also make use of the 
fact that 

+~ 


1 ۲ y 1 1 n 
— —tox® ___ = — سس ہے‎ 
T f e iv| dy = 6. (x9) — à, (a5) xi 


بح سے 


and choose f (xy) = g (مد)‎ = — 17 or f and g in (5.24) as a solution of the 


homogeneous equation, We t arrive at the dcs 


2 
AM z— n p? nn ume u? +4 2]. IE (5.33) 
9x? 
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Thus the symbolic method just discussed, which is based on the analogy with 
the Yukawa equation, is closely connected with approximation method A, while the 
symbolic method based on the analogy with the oscillator equation is closely con- 
nected with approximation method B. As we see, approximation method A reduces 
92 
dC 

'The above also gives the mathematical justification of the assumption as to the 
order of u in method A. With a different choice for the order, we would actually have 


to a series expansion in powers of the operator pet 


2 
to make the expansion of wet ee , which would lead to serious difficulties. 
x 
We wish to thank J. Werle for interesting discussions during the writing of 


this paper. 


Appendix I 


We desire to perform the summation in the expression: 
~ 1 
ue l - I 
dir us سر‎ ru) emos: 1 


REO 
۱ 1 ; 
We represent ô (t) in (A.1) by On / do e“ , Thus: 
7 


1 رہ‎ 21 
epa —io == 
4 E f ev iot ب2‎ 7 ) 7 w?) > 


يم — 


Te 


1 : r 
` Baer IE mr c T res 


-m~ 


1 


q — 

M hoe Ue E چس چ‎ 1 

== 8r2cr fa e-t ch 23 q? — w2 + کرس‎ Te I cos L w2 — q?. 
—q رك‎ : 


Remembering that u = q/c, the above expression can be represented by 


cuc 


1 w? 
Av = do e-iot e-r ۷۳-7 


8z?cr 


—uc T : 
1 : Fist 
jd عم‎ 


Be 
1 
2 gag | dcm a 


c? 
-m . 
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We let: 


-Fuc 
1 1 w? 
OS = —iot Cif en acc z 
2 af نا ہے 1 م67‎ c ES 
Says 


Employing (A.3) and the representation for ‘A given in Table 1.4, we see that 
Av = A+ 3 /4odd (A.4) 


Since A® and 4 are solutions of the equation (L1—g?) G = — bw (x^), then As 
defined above satisfies the homogeneous equation: 


)0- p?) 4 =0 (A.5) 


4^? can be given explicitly. 


+n/2 
1 e incising ch (ur cosp). 


—a[2 


my dro 
An? r 7 
However, since:1 
22 
1 fao ch (a cos g) cos (b sin p) = I, V(a? — 53), 
7 oe 
thus 


e رو‎ Ae 


Ar r 7 


1 This formula is very easily obtained as follows: 
27 


7/2 2n | 
1 : i a cos p + ib sin 9 dg = i d uti ی‎ E 
= dq ch (a cos (و‎ cos (b sin 9) = ہے‎ 6 A pea 

0 0 


| —m[2 
: 27 


ear er یتوم و‎ 
2 3 سا‎ N 3 ein me dg = 3 هوه مه یر‎ 7 (Va? — b$). 
n=0 


n! ^m! 27 n? 


mm=0 
As seen from this proof, the formula holds for any arbitrary complex a and b. 
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Performing the differentiation and distinguishing between the regions inside and 
outside the cone, we therefore have: 


y ) (| ۲ xax? 20 


٦48 — p? 
I, (u 0(| تم‎ xax > ۰ 


- (A.6) 


Substituting this into (A.4) and using A from (2.4a) we find A® explicitly: 


1 Su? 0 qo uo 0 


Iro Ho xa x* < 0 (A.7) 


Appendix II 


From Bessel’s equation and the recurrence relations for Bessel functions we ob- 
tain the relations: 


9 ۱.2۶ (qt) _ . 1 9 ء2‎ )01( 
(Bte) P- -ea (A.8a) 
eA ate Zen g) 


Combining (A.8a) and (A.8b) gives: 


92 | Z, (qt) Z,+1 (qt) 
(5 + 2 oP - تیه و رز + رن‎ 


from which it is easy to show: 


9? igi ZI b ا ریا‎ 
(se) A E ee 2n Mgt due S (A.8d) 


And in particular, 


92 q 
دا‎ 2 (qt) = 3^ 2, (qt) (A.8e) 


Appendix III 


We will examine the series: 


ABIT Zee ES تروع اکب‎ 
1 EY zl ng a] Jm (at) + 


JY A fe 1 z= ] +7 
Tand ےم‎ a an 2,2" ( k | ecc 9 وم‎ 


of this paper). The question of the convergence of the first part‏ 4.24 و و( 
of 4°“ presents no difficulty. From (4.18) and the convergence of Lommel’s‏ 
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: A t 
series, we see that the first part of (A.9) converges to the quantity 11 Ada. 


We let: 


v 1 = ېم‎ 1 c Pus 
P9 3 al توس‎ 2k 2 | §(2l+1-2k) (t 
I d gx à. ° | k L ووو‎ (A30) 


k=0 


Thus, we have: 


٢٢ — oA + m 4d (A.11) 


In order to perform the summations in (A.10), we note from the properties of 
the Euler integral B (و ,م)‎ that: 


oe (i)[ea-m ال‎ 
0 


Inserting (A.12) into (A.10), we have: 


1 
: d ee) ey We دي‎ el ود رک‎ 
deas AER u خی‎ eir. Q (1) = 
i xa وا‎ 7a Jara 29 | اخ کا ام‎ O 
: : 


1=0 k=0 
= 21 : 1 
ارم‎ E 1 21 q? 2? (ys) 7 
ات ود‎ AE E) 6 (t) =  (A.13) 
y 0 
2 : 1 
1 v^ 7 4 22 :] 7 
~ e 113 0 fe 0 T) Ot? +9 0(6). 
= 0 
Reversing the order of summation and integration, we have: 
2 1 
کت ھ2‎ dt J ae — (1 — 7( ا‎ - q? | ۵00۵ (A.14) 
2727 2 ot? 
0 


Putting ۵ (t) in a Fourier representation, we have: 


1 +~ ۱ 
w 1 7 4 
سس‎ fo fiwon Jo E WU nor ce s) لي‎ D) 
E | 
jo 


1 


TÊ 


1 r dt z 
=a | (2 Vr? — 7) ۳ سے‎ Ai 
ہم‎ 0 
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(v = V1—7? has been introduced as the variable of integration). But: 


iv 


1 co 
r : : t 
0 0 


ITA 


Using (A.16) in (A.15) and taking note of the parity, 


co 


i [onn ۰ت‎ ho = 


— Anc? [t| . 
0 
۲ (A.131 
7 EvA) 
1 3 
= کی ہے‎ dv v ول‎ (vt) ړل‎ y2 = q? 


0 


We compare equation (A.17) with the discontinuous Sonine’s Integral (Watson 


1944, p. 415) 
1 o o<a<ß 
J; J, (a Vx? + 22) اھر‎ » [Voi — B3 21 
الا‎ cw دح‎ 
0 a>f>o 


(A.18) 
which holds for any complex z when Req > Rep > — 1. We thus obtain: 


0 o R روک‎ 


: DEL HE pra (A.19) 
ہت ہہ سح‎ ee 2 2 2 42 
A ات‎ TOS و‎ m a) es 


"We will now perform the differentiation indicated in (A.19) taking into account 
the discontinuity of the function being differentiated. Noting that u = q/c, we have, 
finally : 


NE Y E 
۵ (c? t? — r?) + p? پل‎ ۷۷۳ — et) prar) 
47۱1 u Vr? — èt? 


(A.20) 


v 2 I. 1 
T N EE 


t?‏ لع یں 


1 From the relation 


1 ۵ 
foie zl, (z) = 1, (z) 


we obtain 
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Comparison of (A.20) with (4.8) gives: 


A= —2— je. (A.21) 
From (A.11) we thus obtain: 


LEX == E ZI (A.22)‏ 5ت — — Aoda)‏ + 2/۳ —( سر 
using aN "t (2.6).‏ 


Equation (A.22) once again verifies (4.29). At the same time, its derivation shows 
how the series in 1/6 which represents the function A converges. The ‘regular’ part 
of this series converges on the strength of the convergence criterion for Lommel's 
series. The convergence of the singular part-of A —if it makes, any sense to speak 
of it — is of a type similar to the convergence of the series for D and D formed by 
6 functions. This convergence could be checked ا‎ investigating the convergence 


of the series obtained from fa dt’ ‘i d, x 2 TESTA: x— x’) p (t' x) by reversing the order 


یم 


of the summation and integration 5 making use of the properties of the 6 function. 
We may add that the auxiliary function 4, which we د‎ obtained in (A.21) 
1 
2 [t] 
that this function, which differs from zero outside the cone and equals zero inside the 


is related ما‎ ۸۳ by an expression analogous to 4 = اھ‎ easily be seen 


the cone, satisfies the homogeneous equation (CJ — 4?) A = 0. For these reasons, it 


would perphas be more logical to denote A in a different manner, namely, 4 = ۰ 
However, it is not worth attaching particular importance to this, since it is not a rela- 
tivistic function. 


Appendix IV 
Equations (A.8a) and (A.8b) can be written: 


2۶ 259 ٠ا9۵‎ _ EA 
E de imd اب و ۳ مس(‎ TET وت‎ 
1 2 Z (gt) Za (A.23b) 
d OL gy 


If we let Z, (z) = H™ (2), replace t by r, q by iu, and use the definition K, (z) = 


۰ 77۲ 
—iv ۲ 
= S e? H® (iz), these equations may be written: 


2 2 بش‎ p) یگ‎ - 20-D + 1[ l 


3 
3 


+ لست سے بو‎ (A.24a) 


1 9 KD) نے‎ _ , Ku) (A 24b) 


ee 


۶ 2 ad mii 
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from which: 


9? 23 K, (ur) K, +1 (ur) 
شتو هر سس می مار‎ tac rapa 1] یب‎ , A.25a 
E | r or ۸ ) T VAI (v Per ] +1 ( ) 
: 2 ۵ 
We iterate this relation and regard * as the radial variable so that 5 + ma 
— V2. We then conclude that: 
[2(1—1) +1]! Kia (ur) 
(v? —u3)*! Ky (ur) = — ui acc HAE CEN (A.26) 


٦ 
KPATKOE COJIEPDXXAHHE 


JL. 111061875 u T. Ilze6auckuit, O paseumuax 006600001  dyHkuuü, ۴ ۴ 476 
c 206146146۸4 Kaeüna-I'opoona. 


Ians WA ypasmHenus LUreitua-lop;toHa J[Ba aJIBTepHaTHBHPbIe MeTO/IA IIpHOJIEDKeHH51, CBA- 
saHHble c Pa3BHTHAMH B C-l, JIHOO € c. 
TlpoumMeHAA 9TH MeTOJIbl, IIOJIVUeHbI pa3BHTH51 607166 BaxKHbIX oco6enubix dynkinit, TAKHX, 
Kak 4, A, AC). IleiicrBue 270 cBs33aHO C J[ByMS aJIBTepHaTHBHbIMH MeTO/IAMH (DODpMaJIBHOTO 
uurerpHpoBaHusa ypaBHeHus Lureita-loproHa, 62310110۳0۲۲ Ha aHaJIorHH c ypaBHeHHeM YH- 
KaBy, HJIH 7102246 C ypaBHeHHeM rapMOHMUeCKOrO ۰ 

1 7321535 11671 Tpy/laà, 370 HpHTOTOBJIeHMe MaTeMaTHUeECKHX 0060608 هت2‎ IIOCJIe/IOBATEJIb- 


Horo 18011660۹5 Meroga Oüunrreinua-lnudeznsna-loddéwauma B Bonpoce o6 ypaBHeHHAX JIBH- 
)KeHHA B Me3OHOBBIX 60۰ 
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DIRECT INTERACTION AND THE COMPOUND NUCLEUS 
FORMATION IN NUCLEAR REACTIONS 
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The formulas for the reaction cross section are derived taking into account both 
the compound nucleus and the direct interaction mechanisms. The derivation is based on 
Bethe’s old treatment of nuclear reactions. The final formula for the direct interaction reac- 
tion amplitude is obtained in the Born approximation. The compound nucleus reaction 
amplitude is calculated for one level of the compound nucleus. The cross section is built 
up with the compound nucleus, the direct interaction and the interference cross sections. 
The formula is applied to the (n,d) reaction and the results compared with those obtained by 
Thomas. The treatment of this paper takes into account the transparency of the nucleus. 


8 1. Introduction 


Up to the last few years, the hypothesis of the compound nucleus formation 
was the fundamental assumption in the theory of nuclear reactions. According to 
this hypothesis, suggested by N. Bohr (1936) and Breit & Wigner (1936), the com- 
pound nucleus is formed in the first stage of the nuclear reaction. This compound 
nucleus decays in the second independent stage with the emission of the reaction 
products. The basic assumption made is that there is a strong interaction between 
nucleons which results in the mean free path of a nucleon in the nucleus for inter- 
mediate energies )> 50 MeV) being much smaller then the nuclear radius (save for 
the lightest nuclei). For higher energies, this mean free path is comparable with the 
dimensions of the nucleus since the cross section for nucleon — nucleon collisions 
is then smaller. For such high energies the transparent nucleus model of Serber (1947) 
applies. 

However, both theoretical considerations and experimental evidence indicate that for 
lower energies as well, the compound nucleus hypothesis is not as legitimate and as 
certain as previously hoped for. The great success of the shell model based on the 
assumption of nuclear orbits suggests that also for low energies the mean free path of 
a nucleon in the nucleus is probably at least comparable with the dimensions of the 


(249) 
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nucleus. It probably reflects the exclusion principle, the meaning of which was overloo- 
ked in previous discussions (Weisskopf 1950)". The optical model of Feshbach, 
Porter and Weisskopf (1954) for nuclear reactions with neutrons is in agreement 
with this conclusion. ۸ semitransparent nucleus conception for this model was success- 
fully applied to low energy neutrons (0—3 (۰ 


Angular distributions of numerous reactions of the stripping, pickup and knock- 
-out type provide us with direct experimental data which indicate a mechanism for 
nuclear reactions differing from the compound nucleus formation. The main feature 
of these reactions are pronounced maxima of the differential cross sections for small 
angles. All these experiments are to be easily explained by the assumption of a direct 
interaction mechanism, in which the incident particle transmits its momentum di- 
rectly to one nucleon (or a group of nucleons) of the target nucleus. Such an interaction 
for higher energies results from the transparency of the nucleus (e. g., Chew and 
Goldberger 1950). For lower energies (— 10 MeV) the stripping, pickup and knock- 
-out reactions, in the majority of cases, were satisfactorily explained as surface effects 
(Butler 1951, Austern, Butler, and McManus 1953; see also the review article of 
Huby 1953)?. However, the explanation of these reactions probably should not be 
based on the surface mechanism. There are indications that the direct interaction 
involved here is not confined only to the surface of the nucleus, and this fact is in 


agreement with the transparent nucleus model (Dabrowski and Sawicki 1955a, Se- 
love? 19564). 


One has to state however, that save for the above-mentioned stripping, pickup 
and knock-out type experiments, a great number of phenomena, especially reso- 
nance phenomena in nuclear reactions, are satisfactorily explained by the compound 
nucleus hypothesis. So we have two „complementary“ descriptions of nuclear reac- 
tions which can be used for the calculation of the cross section. One is the theory of 
the compound nucleus, the other is the theory of the direct interaction. The mutual 


relation of the two descriptions seems, at present, to be the fundamental problem in 
the theory of nuclear reactions (Weisskopf 1955). 


It is the purpose of this paper to derive the formula for the cross section for nuclear 
reactions, taking into account, simultaneously, the direct interaction and the com- 
pound nucleus formation. Obviously, the reaction amplitude is, in this case, the sum 
of the amplitude representing the direct interaction mechanism and that representing 
the compound nucleus formation. The derivation of this formula is essentially a stan- 
dard one. Such a formula is desirable especially in view of numerous stripping experi- 


1 See, however, the paper of Eden and Francis (1955). 
* Compare also the review article of Dabrowski, Postepy Fizyki, 6, 274 (1955). 
* The author is indebted to Dr. W. Selove for letting him know of his results prior to publication. 


4 See also the results based on the transparent nucleus model for the pickup and knock-out type 
of neutron-induced reactions on Li (Dabrowski and Sawicki 1955). 
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ments recently performed in the vincinity of an isolated resonance of the compound 
nucleus}. 

Recently, R. G. Thomas (1955) has given formulas for the (n, d) and (p, d) 
reaction cross sections taking into account the direct interaction and the compound 
nucleus formation. His elegant derivation is based on the R-matrix theory of Wigner 
and Eisenbud (1947). The R-matrix theory is of course the most correct and most 
general approach in the theory of nuclear reactions. Nevertheless, it is the most 
formal approach. However, the distinction between the direct interaction and the 
compound nucleus mechanism is rather more intuitive. Consequently, Thomas 
utilizes the most intuitive point of the R-matrix theory: the division of the configura- 
tion space into an internal and external region. The consequence of this procedure is 
the treatment of the (n, d) and (p, d) reactions as a surface effect?. 

To get an expression of the direct interaction amplitude taking into account the 
fact that the direct interaction is not only a surface effect, i. e., taking into account 
the transparency of the nucleus, one must base the derivation of the formula for the 
cross section on a more intuitive theory of nuclear reactions than the R-matrix theory. 
The treatment in the present paper is based on the old treatment of nuclear reactions 
of Bethe (1937), in which the direct interaction and all spin consideration are included. 
The derivation is restricted to one isolated level of the compound nucleus. It seems 
that one can take into account several levels without any essential difficulties. In the 
final formula, the direct interaction amplitude is given in the lowest (Born) approxi- 
mation. A better approximation e. g., the distorted — wave Born approximation) 
does not present any essential difficulty, but makes the final formulas more cumber- 
some. The application of the formulas obtained in this paper for tne pickup reaction 
is discussed in 8 3. The results are of course identical with the Born — approxima- 
tion formulas of Thomas. The only difference results from the integration over the 
internal region appearing in the formulas of this paper. 

For the sake of simplicity, the Coulomb effects are neglected in the present 


- derivation. 


8 2. Derivation of the Formula for the Cross Section 


I. Notation and Fundamental Equations 

The following reaction is considered: the target nucleus X, is bombarded by the 
particle P,, the particle P, emerges as the product of the reaction, and the final 
nucleus is X,. We denote this reaction: 


PATS IDE a. D 


1 Berthelot et al. 1954, 1955, Van Pattern et al. 1954, Juric 1955, Jones et al. 1955, Stratton et al. 
1955, Risser et al. 1955, Green et al. 1955, Ward, and Grant, 1955, Rosen and Stewart 1955; Boner et al. 
1955, McEllistrem et al. 1955, Grosskreutz 1955, Wilkinson 1955, Marion et al., 1955, McGruer et al. 1955. 

2 The R-matrix theory is regarded by Thomas as responsible for the compound nucleus mechanism 
only. But any nuclear process occurring in the interior of the nucleus (including direct transitions) can be 
represented by the R-matrix theory (compare Peaslee 1955). 
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We introduce the following notation (all quantities are measured in the center- 
-of-mass system): 
a — channel index (a defines the type and state of the particle P, and the nucleus X لل‎ 
s, u — channel spin (sum of the intrinsic spins: i, — of the particle P,, and ی‎ — 
of the nucleus X,) and its Z component!, 
,ے9‎ — normalized linear combination of the products of the internal wave func- 
tions of P, and X, with the resulting channel spin s, ځا‎ 
l, m — orbital angular momentum and its projection on the Z axis, 


YM,— V. [um] GlumJ M) 9,,, ور‎ ) ۵ 
PY —Y RA, (r,) — wave function of the system in the channel a with the channel 
spin s, orbital angular momentum J, total spin J and its Z component M, 


YP, sc 5 {sl 7 M} VÀ, — wave function of the system in the channel a, 


(jı jg m4 m,ljm) — the vector-addition (Clebsch-Gordan) coefficient’, 

Yin — normalized spherical harmonic’, 

Fa (as Ra) — vector from the c. m. of X, to the c. m. of P, and its polar coordinates, 

RAE, (ra) — the radial part of Wf, (regular and with proper asymptotic behaviour), 

XM, = X, (f, M, zt) — the normalized wave function of the compound nucleus in a 
state with total spin بر‎ its Z component M, and with parity څل‎ 

H = H, + T, + V, — total Hamiltonian of the system, 

H, = Hx, + He, — Hamiltonian of the internal motion of X, and بیط‎ 

Va = Vx,p,— the interaction between X, and P,, 

T, = — (h?/2M,) A, — kinetic energy operator of P, and X, (4, = AF), 

M, — reduced mass of channel a, 


h k, — momentum of P., 
v, — relative velocity in channel a, 
W = E, + W, — total energy, 
W, = Wx, + We, — internal energy of X, and P, 
W,— energy of the compound state, 
E, = h* 12/2 M, — kinetic energy of X, and P, 
The wave function Y of the whole system must satisfy the Schrödinger equation 
(H — W)w — 0. 0 


We assume the following form for the wave function V: 
¥ = 3 Pork رز‎ M3 Cm, XM, . (3) Š 


1 The Z axis is taken in the direction of the velocity of the incident particle. 
2 The behaviour under the time-reversal operation K of all wave functions is: 
K y= )-( ۳ ". 


This representation was suggested by Biedenharn and Rose (1953). 
* Defined as in Condon and Shortley (1935). 
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In Eq. (3), the assumption is made that only one state of the compound nucleus 
plays an important role in the reaction considered. This is well-justified for energies 
near the energy of an isolated level of the compound state (W ~ W). We assume, 
further, that the wave functions Zu, of the compound state are restricted to a finite 
region of space. 

Since 


1 و‎ E W, P Him, == 2 XM, (4) 


by inserting the representation of V/ as given by Eq. (3) into Eq. (2) one obtains for 
C,, the equations: l 


Cu W—-W) = V (asl) < tml ۲۵ - ام‎ Wily, >, 6) 


and for the functions RM,, the equations: 


[(h?/2M,) ddr, + E, — Vasy (ra) - AL )1 + 1(/2 Mari) Rs (ra) Ta = 


EN (e s (ا‎ fdr, dQ, واد‎ — L) VAL, + 


+ Cu dy, [ dr, 49,7, Y, (Vn — LY gas (6) 
where 
L, کے‎ E, Piel, a? - l (7) 
and 
ANIE f dt, 49, Vit V, وی‎ (8) 


The integration f dr, = if; dtx, dtp, extends over the internal coordinates of the nu- 


cles X, and of the particle P,. The integration in the matrix elements > | | > is 
performed over all coordinates. The summation 2” excludes the term {asl} = (a' s' l'}. 
As a consequence of the rotational invariance of V, L, the effective potential ۳ ے‎ 
is independent of M, and the functions integrated in Eqs. (5), (6) have the same 
indices f, M. 
We now devide the wave functions V^, into two parts: YW, representing the com- 
0 


pound nucleus mechanism and V^, representing the direct interaction mechanism. 


D 
The incident wave in the entrance channel a = په‎ is included! in V/,; consequently 
[$5 سے‎ 


the functions V/, contain only the outgoing waves. Thus we have: 
D 
سو وا‎ OP اا بل‎ Gel My vat SE aio 
0 کا‎ FEE 0 
Wow wou LA O) 
0 D 


Jj 
" 


1 This is a question of convenience. We may include the incident wave in ٣ equally well. 
D 
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where the superscripts „inc“ and „out“ denote the incident and the outgoing waves, 
respectively. The division (9) induces, of course, an analogous division of the radial 
functions RE This division is defined by the equations which are satisfied by the 
functions Rs and HOST Namely, we require that 

C D 


[(h3/2M,) 2*۵۳ + E, — Vy — h3 1 (1 + 1)/2M, ral Rasy T, = 


oy NG Voy Va Lota 5 
and 
[(h3/2M,) ddr + E, — ویر‎ - BL ( + 12M, ra] Rely Ta = 
۱ (10D) 


= V (a s) f de, ت80‎ YM, Va — ویک ط‎ 


Eqs. (5), (10C), (10D) are our fundamental equations. All these equations are 
coupled together. As an approximation, however, one can solve the equations (5), 


(10C) and (10D) independently. 
II. Solution of equations (10C) 
We introduce the notation 


V SP VT | (11) 


V, plays the role of the internal interaction responsible for the formation of the com- 
pound nucleus. 


Now, we can rewrite Eqs. (10C) in the form 


[(h3/2M.) آ77‎ Vas 121) + 1)/2M, 77] Rast Ta = 


= (es ny ۶۶۶۶7ٔ7“+ — L) tae (12) 


Let us denote by f and g, respectively, the regular and irregular solutions of the 
homogeneous equations associated with Eqs. (12): 


Jas (r ce 0, 13 
[(°/2M,) 02/07 + E, — Vag — ۸2 ) + 1)/2M, ri] 17 (Ta) (13f) 


۱ Sasl$ 7 2 0, (13g) 
with the following asymptotic behaviour for large values of r: 


; 1 
fag -— sin c ۰. DX lx + à) 7 (14 f) 


1 
زامهق‎ > COS (s E NT lr + bas) ۰ (14g) 
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For large values of r,, for which the right-hand side of Eqs. (12) vanishes (owing 
to the assumed behaviour of ¥4), we have 


M M M 
Rasy fa > Casi زا هل‎ + Dasty Eang" (15) 
For all values of r,, we can write 


Rosy اچ وی‎ Gag Jay d. p Sasl¥ (16) 


where 2a is the solution of the inhomogeneous Eq. (12), which is identical with وروی‎ 


for large values of r,. Obviously, g,, is regular at the origin. On introducing the nota- 


tion 
M M* T 
tee el Yn f asl To» ا پ6‎ YN اوه‎ ۸ (17) 
we may write 
M ھ‎ pM M (M 
DD = Gast 80 PETI TE Pasig C ; (18) 


In order to determinate hers we multiply Eq. (13f) by Agee Ta» subtract it 
6 


from Eq. (12) multiplied by fazy» and integrate the result over r, from 0 to oo. Using 
the value k, of the Wronskian W ج)‎ f), we get finally 


M = — (cul Ava) Vas OI (19) 
where l 
Yast = (2/Vho) < Fas, | Val Zm > (20) 


Note that y,,; is independent of M (a consequence of the rotational invariance of 


V.) and is a real number (a consequence of the time-reversal symmetry of V.) (compare 


Biedenharn and Rose 1953). 
We now detrmine the coefficients anny: Let us write 


ibe p" Any (21)‏ = وه 


From Eqs. (9) we see that 
سا کا‎ pee ہے‎ yout 
c 


C 


contain the outgoing waves only. On using the well-known expansion (e. g. Blatt 


and Biedenharn 1952) 
D Em 0 eon d) و‎ = i (ko Ty 3 Uo JIM} (Solo Mo 0 | jM) (21g + 1۳4 
3 1 : 
Lr M و‎ d (er — 31) t el i (kro — 2) | ما‎ (22) 


one easily obtains 


AM, = 227 کې(‎ easi (sy 1 pty 0 | JM) (20 + 1( 8,,0,,. (23) 
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Eqs. (18) and (21) vield finally 


)24( ما اف ورو 3M) YM,‏ 8# سس 
with 84, end AM, detemined by Eq. (19) and Eq. (23), respectively.‏ 


III. Solution of equations (10D) 


The calculation of this part is restricted to the lowest Born approximation. We 
introduce the additional assumption 


f dx, dQ, YM, (V, — L) VM, = 0 for! ¥ بر[‎ s + 53 (25) 


In the most important case of spin-orbit coupling, this assumption is exactly 
the case. Generally, one can treat this assumption as the neglect of non-diagonal matrix 
elements of V, in comparison with the diagonal ones. Owing to assumption (25), 
in the Born approximation the direct interaction contributes only to the reactions 
(a + aj). In obtaining the reaction amplitude for the real reactions (a + ag), assump- 
tion (25) is unnecessary. 


Using Eq. (25), we may rewrite Eq. (10D) in the form: 
[(^3/2M,) 2/0۳ + E,— h? 1 (1 + 1) /2M, r?] Rosy ج‎ 
Vag Ry ١)۰ (26) 
In the Born approximation, we put گا = لا‎ on the right-hand side of Eq. (26) 
(note that Ps contains the outgoing waves only, consequently we set Ras asiy = 0 on the 
right- ew an of Eq. (26)). We thus obtain 
[(^*/2M,) 0*0۳, + E, — h? 1) + 12M, 7] RM, rq = 
= (1—4,,) f FON YM, (V, — L.) pin, (27) 
The procedure leading to the solution of Eq. (27) is analogous to that used for Eq. 


(12). The functions f, g corresponding to the functions کر‎ g of section II are 
DD 


fu (ra) = ka Taji (ka Ta); Eal (ra) بی یج‎ ka Ta n (ka rj); (28) 
D 

where Jj» ^, are the spherical Bessel functions (Schiff 1949). The quantity analogous 

to A in Eq. (21) vanishes, since W, contains only the outgoing waves. As the result, 


D 
we get the فان‎ asymptotic form of سه‎ 


Ray > Tur ٦ 2 | Wine > را ا‎ A n) )1- ۵ a0: (29) 
ae A(P is the spherical Hankel function of first kind (Schiff 1949). 


1 The part of this equation with L, is satisfied automatically. | 


—— 
B 
E s 
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On using Eq. (29) with the help of the well-known expansion of e ara and the ortho- 
gonality of the Clebsch-Gordan coefficients, one easily obtains the asymptotic form 


of V: 


D 
/ z^ 2i {asu} (vol v.) * eara نو‎ d AP (a Su; 050/0(« (30) 
D 
where the reaction amplitude 
ae (asu, دس سی‎ be 
— Cah’) (M, Mokate)” دب‎ Pay [ V, | در‎ (1 baa). )وا‎ 
It is well known that, in the matrix element of Eq. (31), one may replace V, by Va 
(e. g., Schiff 1949). 


IV. Determination of Cy, 
To determine C,,, one must insert into Eq. 9 the functions 


M 
شه ا6 7 = وزیا‎ 


The function D, is given by Eq. (18). A similar expression is obtained for the function 
C 

HA On applying the Born approximation, we put (similarly as in section IIT) on 

the right -hand side of Eq. (5) 


M uy M 
Pag, — Pe. IE - 0. 

This procedure is based on the fact that in the vincinity of the resonance the compound 

nucleus mechanism probably plays the most important role in the reaction. On taking 


into account the function ¥ in Eq. (5), one would get an higher order process, namely 
D 

a direct transition followed by the formation of the compound nucleus, which finally 

decays by emission of the final reaction product. Such processes are neglected in our 

approximation!. 


On introducing the expression (18) for کی‎ into Eq. (5) one gets 
Jc و‎ 


1 (hw) Ass, Yes - Wa + T2), (33) 


h ڑاسے‎ —Sol 
where 


r= V {asl} yis وک‎ (33) 


and where W, now denotes the former value of W, from which the real quantity 


1 xc 5 1 
6 = 51 2 {as lj Yast 2(h Va) 7% > xm | Va — Lal Gay, T Se 


has been subtracted. ^ 


1 If we were to take into account y in Eq. (5), a more precise treatment of Eq. (26) would be requi- 
red. D 


af 
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V. Cross Section 


By inserting the expression (32) for Cy into Eq. (19) one gets 


Jets 
M c a ky oo) X} a3 (s وا‎ fy 0 | JM) + 6 
امعس | وا‎ 
Yad Yann V ا (2 کل‎ Oy, (35) 


Using this expression for 8%, in Eq. (24) we get finally 
rove SNA sut (ulu) هی‎ r -3 Dosu 4° (aspi, مه‎ So Lo), (36) 
C 
where 
A (asit, agsoto) = in^ ko رز‎ ly mM] Ql, + 1)^ (slum | JM) ۰ 
(soloko 0 | JM) [0,,,0,, 0, - یه رک‎ I (Q,). (37) 


The scattering matrix S is 


Shia, = تا‎ ۲ du, On, gaa 7 eit د‎ )38( 


This is the well-known expression for the scattering matrix for one level of the compo- 
und nucleus (see, e.g., Wigner and Eisenbud 1947, Blatt and Biedenharn 1952). 
On joining together Eqs. (30) and (36) we get finally 


Day, ٥٦ (39)‏ بر دوو 

where the reaction amplitude 
مور کٹ‎ (40) 

The cross section for the reaction ولاوکمه‎ > asu is 
do (asp, AySoM)/dQ, = [do(asu, ہ۸(/4۵لمتمہ‎ compound + 
+ [do (asu, مه‎ So ug) dO] pis. + [do(asp, AgSolto) | dal interference, (41) 
where 

[do (asp, aos) dO Jcompouna = | 4° (aspi, رڈ|(وطلوتوہ‎ (42) 
[do (asp, soto) dR d] Direct = | AP (asu, a95oHo) |? (43) 


[do (asų, a959/40)/ dO. tnterference کت‎ i (asy, I) AP (asu, A) + cc. (44) 


From ihese expressions, one gets the cross sections for the reactions Qg$g — مه‎ 
and a — په‎ by summing over the final channel spin states and averaging over the initial 
ones. These calculations for the compound nucleus cross section were performed by 
Blatt and Biedenharn (1952). The corresponding calculations for the direct interaction 
and the interference cross sections may be performed for concrete reactions. 
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8 3. Discussion 


At present, the most important case in which the formulas obtained in 82 may 
be applied, is the case of the stripping reactions (d,p) and (dn), and the inverse pickup 
reactions (n,d), (p.d)! for energies near the resonance energy of an isolated level of the 
compound nucleus?. We can discuss, for example, the pickup reaction (n,d). In the 
direct interaction amplitude we take into account the interaction between the incident 
neutron and the picked up proton only (we replace V, by V, = V,,in Eq. (31))?). 
We assume this interaction in the zero-range approximation (Horowitz and Messiah 
1953). 

It is possible to perform the radial integration appearing in A”, if we assume that 
the proton (which is picked up in the reaction under consideration) moves around the 
final nucleus* under the influence of the effective potential Vj,1, (r,) (J, is the sum of 
the spins of the final nucleus and the proton, 1, is the orbital angular momentum of 
the proton bound in the target nucleus, r, is the distance between the c.m. of the final 
nucleus and the proton). Let us denote by Rj,1, (r,) the radial part of the wave function 
of the proton in the state j,, lp- We use the result of Daitch and French (1952): 


1 dr, 15 Jip (yrs) Riply (Tp) = 


(ur) Ry, (7), 09)‏ سر ور لو مل سه 
i 0‏ 
where a is defined by the condition‏ 
Vjy, (r,) = O forr, >a, | (46)‏ 
and?‏ 
Go‏ ذرامنظ Viplp me f dr, r5 Jip Riplp Viply | Ji dr, rg Jip‏ 
ô 0‏ 


the remaining quantities have the same meaning as those introduced by R. G. Thomas 


(1955). 


1 In our treatment, the theory of all four of these reactions is identical owing to the neglect of 
Coulomb effects and owing to the ,,reciprocity theorem for nuclear reactions“ (Blatt and Weisskopf 1952). 

2 See the experimental works quoted in footnote’. 

3 Thus we assume that the compound nucleus formation is the most important mechanism of inte- 
raction between the neutron and the final nucleus, this being in agreement with the conclusion of Francis 
and Watson (1954). ۳ ۱ 

4 This assumption is supported by the arguments of Peaslee (1955), who suggests that the direct 
interaction mechanism is of importance in cases with similar initial and final configurations of the nucleus. 

5 In the simplest case of the rectangular potential well Viply(rp) = — V°iply and radius a, we have 


ape,‏ د 
Viyly = Pipl:‏ 
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Making use of Eq. (45) one can easily calculate 4 ! as well as the direct interaction 
and the interference cross sections. We give here neither the calculations, nor their 
results, since the calculations are rather standard, and the results are, of course, identi- 
cal with the Born approximation formulas, of R. G. Thomas, multiplied by the factor 
[1 — ^? (MV;) (k3, +a] and [1 — h? (M7; p) ! (2, + a?]? in the case of AP, 
(dol dO.) snrerference and 0 respectively. This factor manifests the tran- 
sparency of the nucleus in the case for which a one particle model of the target nucleus 
is assumed. This factor causes one of the zeros of the Butler angular distribution 
to disappear (Daitch and French 1952). In the case of the ?Be (p,d) *Be reaction, one 
can explain the experimental angular distribution of deuterons only if one takes into 
account this factor (Dabrowski and Sawicki 1955a; Selove 1956). 

It is interesting to note that the interference cross section enables in principle 
the signs of Ya to be determined. 

In the present paper, only the first (Born) approximation was presented. One can 
go further, of course, on applying the Born distorted-waves method in a manner similar 
to that used by Thomas or by using the modified Born method given by Tobocman 
(1954). This further approximation seems to be important because of considerable 
distortion of a plane wave for energies near the resonance of the compound nucleus. 
The next stage is to take into account a number of resonance levels of the compound 
nucleus. 

It is evident that the formulas of the present treatment are suitable for calculation 
of the spin polarization of deuterons, and of the angular correlation of radiation 
emitted by the final nucleus after the pickup reaction. 

Using methods similar to those used for the stripping and pickup reactions, other 
reactions, e. g., (n, n), (n, p), (n,a) may be investigated with the help of the formulas 
obtained in $2. It seems that the formulas of the present paper may be especially 
useful for higher energies because of distinct transparency of nuclei for such energies. 
Of course, for energies that are, so high that the different compound nucleus levels 
coincide with each other, the interference cross section vanishes in the case where 
the signs of y,, are distributed at random}. 

I am grateful to Mr. J. Sawicki for interesting discussions on this topic. 


KPATKOE COJIEPXKAHHE 


Bosdeticmsue npa wo u 052030601046 6760 AOpa e nOepuwx peaKYyUAX. 

B ۳26016 مته 6ه صتعھ‎ dopMysrmr nomepeunHoro ceuenusa DIA SJIepHbBIX peakIIHii, 55 
BO BHHMaHHC KAK MeXaHH3M CJIOXKHOTO Apa, TAK H BO3JIelicrBHue 7025۸0. ر0صاطظ‎ 3707 OCHOBBI- 
BaeTCH Ha CTapHHHOM 010106066 Bethe’a orHomreuua K "HepHBIM PeaKIHAM. OJakJIOUHTeJIbHas 
QOpMyJlà 3MIUIHTy/bI peaKIHH BO3JIeiicTBHeM Ip3MO nojryueHa B mnpu6]ikeuuu BopHa. Am- 
11712115718 peaKIMH IIOCDÉZICTBOM CJIOXKHOTO SUIDa HCUHCJIeHa JIA OAHOTO ypOBHA 0 
aapa. Ilomepeunoe ceuenue ABNAeTCA cy. 


MMO 11011616111010 06۷2 [3 62141117 17۲ 6 
11021411010 Apa, 1101161061187010 ceueHHs 52031167۶ 11۳17110 a TAKKE 5776 HHTep- 


PepeuumoHHoro ceueHus. POpMyJIbI HpHMEeHeHBI K peaknuu (n, d), a ۳6311811 CBepeHbI 


c peayzbraramH Thomas'a. تسج هه‎ Hacrosmieit pa6orbi YUMTHIBaIOT IpOSpauHOCTb 4 


A. TomOposcrnun, 


1 For further discussion see the paper b Thomas (1955). 
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A NOTE ON THE DIFFUSION CONTROLLED PROCESSES 
IN LUMINESCENT SOLUTIONS 


By A. JABEONSKI 


Physics Department, Nicholas Copernicus University, Torun 
( Received March 3, 1956) 


A physical meaning is given to a constant 4,;, introduced by the writer in his earlier 
papers. 


The aim of this paper is to supplement two earlier papers of the author on lumi- 
nescent solutions (Jablonski 1954, 1955). These papers deal with some processes such 
as the quenching of photoluminescence or self-depolarization of photoluminescence, 
due to the interaction of excited molecules with some of the adjacent molecules; we 
shall henceforth call the latter “disturbers“. 

In the first of the above papers a luminescent solution is considered as a system 
consisting of excited luminescent centres! immersed in a solvent. All the centres are 
divided into groups according to the manner in which their shells are occupied by 
the distrubers (in our original paper we used "quenchers^), a group of centres with no 
disturbers, a group with a single disturber in the first shell, and so on. If a disturber 
in a particular centre belonging to a given group moves from one shell to another, the 
centre will then belong to a different group, and we may call this change of group 
a "transition of a centre from one group to another. l 

The probability of such transition is assumed to be proportional to time. Thus, 
the probability that a transition from group i to group j takes place during dt can be 
written as پ۸‎ dt, where A; (the transition probability) does not depend on time t. 
There is no doubt that 4;; must be connected in some way with the diffusion coefficients 
of luminescent molecules and disturbers in the given solution (Jabłoński 1954). An 
attempt to establish such a relation is made below. 


1 By the ,,luminescent centre“ we mean a luminescent molecule together with the molecules sur- 
rouding it. Different models of a luminescent centre were put forward by the present writer. The 
simplest one consists of a luminescent molecule with an active sphere. The more refined model is 
conceived as a luminescent molecule envelopped by „shells“ formed of monomolecular layers of the 
solvent. The shells may contain some disturbers. To be sure, even this refined model is merely a 
crude approximation to reality. 
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Let us consider, first, the case of the simplest model of a luminescent centre, 
viz. a model consisting of a luminescent molecule with an active sphere, in which there 
may be, at most, one disturber present. In this case, all the centres have to be divided 
into only two groups: the centres having one disturber within the effective sphere for- 
ming, say, group 1, and those with no disturbers, belonging to group 2. Let À, denote 
the probability of a transition of a centre from the group 1 to group 2, and وه‎ that 
of the reverse transition. Clearly, ی۸‎ dt is the probability that a disturber enters into 
the active sphere of a given centre from outside during dt. The process of diffusion of 
outer disturbers into the active sphere may be described in a way identical with that 
given by Smoluchowski (1917) for the process of coalescence occurring in colloidal 
solutions. Let R be the radius of the active sphere, and D the sum of the diffusion coeffi- 
cients of the luminescent molecule and the disturber. Provided that the condition? 


2RYc 


VD > ۲ (t is the mean life of the excited state of the luminescent molecule) 
7 


is fulfilled, Smoluchowski's expression for the probability of an encouter per unit 
time can be used, namely: 


W = 4x DRN, (1) 
N denoting in our case, the concentration of disturbers. Obviously, W = ,ړڅ‎ and thus 
A, = 4r DRN. (2) 
Hence (cf. Jablonski, 1954) 
۸ Aa 
Ay ہچ گے ہے‎ No A (3) 


where v is the volume of the active sphere, and v = Nov. 

The dependence of D on temperature T and viscosity رہ‎ may be given approxi- 
mately if the diffusion of luminescent molecules and that of disturbers is treated as 
diffusion of spherical particles with radii رہ‎ and r, respectively: 
eTA el 1 
` Orn Sins j (4) 
The expression (2) can thus be written as 


0 3r 


)5( چا 1 3 


2 According to Smoluchowski (1917) the probability V, that a particle situated anywhere outside 
of the active sphere in the volume ۲ reaches the sphere before t is given by: 


The expresion (1) is valid only if the second term in the bracket of the above equation may be dropped 
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or, remembering that N = کی‎ 3» 
v 4zR3' 
i 9۹۷۹ہپ‎ "۹ 1 6 
? 3 1 perl ۷ 


In the case of the shell model, the situation is somewhat more complicated. Let 
us assume that the radii of the consecutive shells differ by A R. The volume of the 
n-th shell is then approximately 4x R? AR. Of all possible transitions, the most pro- 
bable are those in which one disturber moves from one shell into the neighbouring 
one, say, from shell n into shell n + 1. Transitions with | An | > 1 are much less 
probable, and, presumably, may be neglected. The probability that a disturber present 
in shell م‎ + 1 enters shell n may be considered as approximately equal to the proba- 
bility that a disturber from ,,outer space“ reaches the sphere with the radius R, by 
diffusion, the concentration of disturbers in „outer space“ being equal to the ,,con- 
centration“ of disturbers in shell n + 1. If there are J disturbers in shell n + 1, this 
concentration must be assumed to be // )4 R? AR). 

Let us now consider one of those groups in which the centres have / disturbers in 
the shell n + 1, and & disturbers in shell n, and second group differing from the 

first in that the centres have k + 1 disturbers in shell n and J — 1 disturbers in shell 
n+l. 2e 

Let us denote by å (k,l > k + 1,1 — 1) the transition probability of such 
a transition, in which the positions of all but one of the disturbers of a centre (indicated 
in the brackets) remain unchanged. By taking (2) into account, we have 


4zDR,l DI 
Mibi: Rak dod dia per TD (7) 


. In the state of statistical equilibrium the following condition must be fulfilled 
(Jabłoński 1954): 


WEL E و‎ p مو‎ met (o, —1 4 Dp, E 25) 


where p, (k, l) and p, ) + 1, 1 — 1) are the relative abundance of the corresponding 
centres in the state of equilibrium of the system (solution). This abundance is given 
by the product of the following two probabilities: 1) probability that there are k 

' disturbers in shell n; 2) probability that there are J disturbers in shell n + 1. If we 
assume that the above probabilities are given by Poisson distributions (cf. Smolu- 
chowski, 1917), we have 


NS و‎ NH of vi 
Pe (k, D) = e NOn tn +1) m dis (9) 


where v, and v, ی‎ are the volumes of the corresponding shells, and N, the concen- 
tration of the disturbers in the solution. 
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Similarily : 
EI gene: D 
1p) m OR CiU SAD we 10 
pe ) + 1,1 1) e ې م) ما کے‎ 1(!)1-1(1 ( ) 


From (8), (9) and (10) we obtain 


Alo k-LI-1 PELID _ j ay 
AK ply 0-0 کے‎ k+1 ۱ 


Since 


2 -2 
ےت‎ ; 
Un+1 Rari Rn ۱ 


eq. (11) can also be written as 


2 — وب‎ 
Aik اا لل‎ (RO 1 او اسر‎ E 
:ا)7‎ + 11-1 >k, l) Rau] ktl kl R, 
1 By means of Eqs. (7), (11), and (12), the transition probabilities can be calculated 
approximately for the shell model of a luminescent centre. - 4 
4 
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LETTERS TO THE EDITOR 


EE عم وم د د دځ‎ 
ON SURFACE TENSION OF IDEAL SOLUTIONS 
By Jan ٦8ء‎ 
Department of physical chemistry, University of Warsaw 


(Received April 10, 1956) 


In connection with Puchalik’s work (1954) a letter to the editor was published by 
Mierzecki (1955) in this journal. An equation was presented, relating the surface 
tension of the solution to mole fractions of the components. The derivation was based on 
the assumption that the parachor of the mixture is additive to the parachors of pure com- 
ponents, and that the molar volume of mixing is equal to zero. That equation was 
suggested for determining the surface tension of an ideal solution, and it was claimed 
that the logarythmic Szyszkowski’s equation (1908) is not adequate for this 
purpose. 

It has to be emphasized that there is no theoretical background for the opinion 
that the parachor of an ideal mixture is additive, neither any indication is known in 
` what well-defined class of solutions it should be so. The appropriate equations for the 
surface tension of an ideal solution were derived by Szuchowicki (1943) and by Belton 
and Evans (1945), and verified experimentally. The derivation is quoted by Gug- 
genheim in his book on liquid mixtures (1952). It may be shown by methods of 
statistical mechanics, that the additivity takes place, but for the exponents of quan- 
‘tities which are proportional to surface tensions 


0 0 
)1( اج werp‏ + ای | مهم - 1 -[ exp‏ 


Here, y is the surface tension of the solution, y? that of pure component ”i” (i = 1,2) 
and a is the area of surface per molecule; x stands for the mole fraction of the com- 
ponent "2". Guggenheim writes: "this is the simplest and most symmetrical rela- 
tion obtainable for the surface tension of and ideal mixture”. Making the substitu- 
tion 
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(yi — y) a 


= Á 1 2 
T 9 (2) 


exp 
we recover Szyszkowski’s equation in the following form 
kT 
y = yi — — log (1 3 43): (3) 
a 


Szyszkowski’s logarythmic equation is therefore directly applicable to ideal mixtures, 
provided that the constants may be interpreted in terms of surface tensions of pure 
components and of the quantity kT/a, as required by equations (2) and (3). 

The amount of work done on the surface tension of solutions must be taken into 
consideration. Not only ideal mixtures, but even regular solutions have been discussed 
and the respective equations have been derived by Prigogine (1948), Prigogine and 
Defay (1950), as quoted by Guggenheim (1952). In addition, the exact statistical 
theory, expressing the surface tension in terms of molecular distribution functions, 
was formulated by Kirkwood and Buff (1949). 

Even the concept of parachor has been severely criticized (see, e. g., Partington 
1952). Although in the case of a pure substance some approximate justification for 
McLeod's rule have been found by Fowler (1937), no similar considerations are 
known in the case of solutions. It seems that the use of the parachor offers no advan- 
tage in the presentation of experimental data on concentration dependence of surface 
tension at constant temperature. The introduction of an arbitrary exponent, ما‎ as 
Puchalik has done, is not convincing, because the value of k is chosen arbitrarily, 
in order to obtain the best possible agreement with experimental data. Besides, it 
should be possible to represent any concentration dependence of any equilibrium 
property of binary mixture by the equation used by Puchalik. It is quite obvious that 
Mierzecki is right when he states that the exact additivity of surface tensions cannot 
take place if parachors are strictly additive. 

I agree that this problem, as it concerns the influence of association, is as yet 
far from being solved. But much cannot be said about the associated solutions on the 
basis of knowledge of the behaviour of ideal, or even regular, solutions, because 
these must be formed by normal liquids as components. I can hardly do better that 
quoting the opinion of H. S. Green (1952), who writes: "There is an abundance of 
evidence that surface phenomena are very different in polar and non-polar fluids. 
In non-polar fluids the depth of surface zone is probably never much more than a few ` 
molecular diameters. In polar fluids, however, the depth of the surface zone is often 
great enough to be observed by macroscopical methods. The experimental evidence 
and probable explanation has been reviewed by Henniker (1949)...” 

Finally, it is probable that some indications, concerning at least the influence 
of the existence of associated complexes, may be found on the basis of a model similar 


to that called by Prigogine and his coworkers (1954) the ideal associated solution 
model. 
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PbSe AND PbTe INFRA RED DETECTORS 


By Tomasz 61 


Institute of Physics of the Polish Academy of Sciences, Warsaw 
(Received April 18, 1956) 


As far as could be ascertained the PbSe photo-conductive and photo-voltaic layers 
processed by the evaporation method loose their photoelectric properties when exposed 
to the atmospheric air. 

The author developed the method, which renders the PbSe layer the property 
of maintaining its photo-conductive and photo-voltaic sensitivities in the atmospheric 
air, which is obviously a very valuable characteristic. 

The PbTe photo-sensitive layers have also been produced by the similar method. 


It had been the earlier study! of photo-sensitive lead-sulphide layers by Sosnowski 
and others (Sosnowski et al. 1947, 1949, Starkiewicz et al. 1946) that has served 
as a starting point and a guide in developing, in our laboratory, of photo-sensitive 
layers of other substances (TlS) (Wolska 1953), PbSe (Checiüska 1953), PbTe 
(Checinska et al. 1954), and improving those of lead-sulphide type. 

It was the success in developing the so called „open“ PbS type of a photo-sensi- 
tive layer, which, contrary to the earlier one, retains indefinitely long the photo- 
electric sensitivity in the atmospheric air (Piwkowski 1953), that led the author to 
` believe that the same property could be rendered to the layers processed of other 
semi-conductive materials (1953). 

During the last few years this belief has been justified in case of T1,5 (unpublish- 
ed — the author) and PbTe (1954). 

As regards the processing methods of PbTe and PbSe photo-conductive layers, 
information is scarse. Some details of manufacture are given by many authors 
Simpson (1947), Treu; Checinska 1953, 1954, Moss 1953, Moss et al. 1948, Kren- 
tzien; Bode et al. 1954, Blackwell et al. 1948, Checiüska et al. 1954, Simpson et 
al. 1951, Simpson et al. 1948). 

As far as could be ascertained, the PbSe layers processed by the above methods 
loose their photo-conductive and photo-voltaic properties in thé atmospheric air, 
whereas PbTe layers, as has already been stated (1954), remain sensitive. 


1 All that is presented in this paper is concerned exclusively with the evaporation method. The 
chemical wet deposition method is not considered. 
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Quite recently, the author of this paper succeeded in developing PbSe „open“ 
type of a layer, whose sensitivity both photo-conductive and photo-voltaic are not 
affected by the atmospheric air, a feature which is an obvious advantage of this layer 
over one that has to be enclosed in a vacuum tight envelope, as it furnishes additional 


possibilities for experimentation and application. 


Processing of the photosensitive „open“ layers 


The method of processing is quite simple. About 2 milligrams of lead-selenide 
containing 5 percent of lead over stoichiometric proportion is introduced into a glass 
envelope. After evacuating it to high vacuum, the substance is condensed on to the 
inside surface of the wall opposite to that bearing graphite electrodes, and then reva- 
porated on to the latter surface. 

The layer obtained is then baked during a few minute time in the atmospheric 
air at 300—400 degrees centigrade, and then cooled down to the room temperature. 
The layer so obtained shows photo-conductive and photo- voltaic properties in vacuo 
and in the atmospheric air. 

It is of interest to note that though a certain number of sensitive layers were 
obtained in condition of high purity, with the liquid air trap between the envelope 
and the diffusion pump on, a few others, just as sensitive, were prepared without the 
latter precaution. 

It appears that photo-voltaic effect of the layer can be enhanced by revaporizing 
the substance in such a way as to make the particles travel in beams inclined to the 
surface provided with the graphite electrodes. The photo-conductivity however, as 
far as could have been ascertained, does not depend on such condition. - 

The procedure employed by the author is the following one. Firstly, a sensitive 
layer is obtained in a hard glass Dewar type envelope. Then, the part of the wall 
which is opposite the layer is removed. In its place a KBr window is fixed with the 
Apiezon lac, to the ream of the opening. The envelope is then evacuated and sealed 
off. The PbTe layer possessing both the photo-conductive and the photo-voltaic pro- 
perty, at low, temperatures, can also be produced by this method. Some details re- 
sembling those described above were used in manufacture of lead-sulphide layers by: 


Berlag, Rumsz and Strachow (1955), Gérlich (1947), and Schwabe (1955). 


Characteristics of the „open“ layers 


Though sensitive layers can be obtained for various widths and separations of 
the colloidal graphite electrodes, the most advantageous for our spectro-photometer 
optical system are electrodes 1 millimeter wide, the tips of which are separated by 
a fraction of a millimeter. 

Variation of resistance of PbSe layers in ohms with temperature in the dark and 
when illuminated by a constant intensity of radiation, of 0.08 Watts/cm. sq., proceed- 
ing from an incadescent, lamps is given in the following table ۱ 


PbSe and PbTe infra red Detectors DIB 


Room temp. | Sotia. CO, temp. Liquid air temp. 
Dark. Illum. Dark. | Illum. Dark. | Illum. 
PbSe Nr. 5 0.45 x 108 ۱ 0.417 x 108 7 — | سے‎ 25 x 108 | 16.8 x 68 
PbSe Nr. 6 0.15 x 106 0.14 x 108 2X 108 1X106 5X105 | 0.21 x 108 
PbSe Nr. 7 0.34 x 10$ 0.31 x 105 —— سے‎ 60 x 108 UDR IUS 


Variation of photo — E. M. F. with temperature, in millivolts. 


Room. temp. Solid CO, temp. Liquid air temp. 
43 z 14.0 
PbSe Nr 6. 12.0 | 112.0 214.0 
PbSe Nr. 7 125.0 | r 937.0 


The overall reactivity of the PbSe Nr. 6 cell, at liquid air temperature (700 1ممع‎ 
chopping frequency) amounts to 1204 mirco volts/micro watt/mm. sq. and the smal- 
lest detectable intensity of radiation to 20 x 10-170 Watts/mm. sq., consisting of 
0-7 microns portion of the spectrum of radiation, proceeding from a black body 
at 220 degrees centigrade, and incident on the modulating disk. 

For the PbTe layer the corresponding figures for the photovoltaic effect were: 
52000 micro volts/micro watt/mm. sq. and 6.26 x 1011 watts/mm. sq. 

No shield was provided that would eliminate radiation proceeding from the 
surrounding objects at room temperature. 

The wave-length limit of the spectral distribution of the photo-conductive and 
photo-voltaic responsitivities, at liquid air temperature of the layers, were; for PbSe 
layer (KBr window) and for PbTe layer (quartz window), 6.5 microns and 4.5 microns 
respectively. 

The author is greateful to Professor L. Sosnowski for his interest shown in the 
course of this work. Wawe-lengths limits of responsitivity were determined jointly by 
I. Filitski and the author. Thanks are also due to Mr. Wardzynski for providing the 
author with the KBr window material. 


KPATKOE COJIEPPKAHHE 


T. IIunkosBcxkuit, J[emexmopu uudpakpacuux ayueŭ uz PbSe u 6 

Hackombko MOHO 6۶70 y6equTca doro-npoBo;pmuue H doro-Bo/nbTaHuecKHMe COM H3 
PbSe nomyuaembie 115701011 BBIIapHBaHHs2 B BaKyMe 116۵2۲۵1 (oTo-uyBCBHTEJIBHOCTb B ATMO- 
CcoepHoM 6 : 

Paspa6oraHo ہہ‎ aKTHBaIAÎ TOHKUX CJIoeB H3 ,PbSe, koropbre co;repoxxunaror doro-uyBcr- 
BHTeJIbHOCTb B aroMrocepHoM 203876. OYeBHIHO YTO 9TO 080116180 OYCHb ۰ 

Tlogo6upm o6pasoM nojyueHo cou 3 ۰ = 
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